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ABSTRACT 

We  compile  the  most  recently  measured  and  calculated  cross  sections  for  pro¬ 
cesses  relevant  to  annihilation  of  positrons  in  solar  flares  and  present  the  final 
evaluated  cross  sections.  We  calculate  the  511  keV  annihilation  line  spectrum 
and  relative  strength  of  the  3y  continuum  for  each  of  the  annihilation  processes, 
using  a  Monte  Carlo  code  that  simulates  the  thermalization  of  and  positronium 
production  by  ~MeV  positrons.  We  calculate  the  thermally-averaged  annihila¬ 
tion  and  positronium  production  rates  and  the  positronium  quenching  rates.  We 
apply  the  results  to  four  specific  environments  (fully  ionized,  neutral,  partially 
ionized  and  a  non-LTE  model  of  the  quiet  solar  atmosphere),  calculating  the  rela¬ 
tive  strengths  of  each  process  and  the  combined  total  annihilation  line  spectrum. 
The  results  are  compared  with  data  obtained  recently  from  the  high  spectral- 
resolution  detectors  of  RHESSI.  We  find  that  positron  annihilation  in  solar  flares 
can  occur  in  a  wide  variety  of  environments.  The  cross  sections  presented  here 
are  also  useful  for  evaluation  of  positron  annihilation  in  other  environments  such 
as  the  interstellar  medium. 

Subject  headings:  Sun:  Flares  —  Sun:  X-rays,  Gamma  Rays  —  Atomic  Processes 
—  Line:  Formation 


1.  Introduction 


Positron  annihilation  radiation  from  solar  flares  was  first  observed  by  Chupp  et  al. 
(1973).  Subsequently,  more  detailed  observations  of  the  line  from  many  solar  flares  were 
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made  with  the  medium  energy-resolution  Solar  Maximum  Mission  (. SMM )  Nal  spectrometer 
(Share  et  al.  (1983,  1996);  see  also  the  flare  catalog  of  Vestrand  et  al.  (1999)).  The  early 
observational  data  included  the  total  yield  and  time  history  of  the  radiation  (Share  et  ah 
1983;  Murphy  et  ah  1984)).  Recently,  the  Reuven  Ramaty  High  Energy  Solar  Spectroscopic 
Imager  ( RHESSI)  observed  the  solar  annihilation  line  for  the  first  time  with  high  energy- 
resolution  Ge  detectors  (Share  et  al.  2003)  allowing  detailed  investigation  of  its  spectral 
shape. 

The  line  has  now  been  observed  by  RHESSI  in  several  flares  (Share  et  al.  2004)  and 
has  been  found  to  be  broad  (4-8  keV),  consistent  with  annihilations  in  an  ambient  ionized 
medium  at  temperatures  above  105  K.  Except  possibly  for  the  first  4  minutes  of  one  flare, 
the  measured  continuum  from  positronium  was  unobservable,  indicating  that  the  density  at 
which  most  annihilations  occurred  was  greater  than  1014  H  cm  3.  The  width  of  the  line 
narrowed  in  2  minutes  to  <1  keV  late  in  this  flare,  consistent  with  annihilation  in  ionized  H 
at  T  <  104K  and  density  >1015  cm-3.  These  observations  suggest  a  highly  dynamic  flaring 
atmosphere  at  chromospheric  densities  that  can  reach  transition-region  temperatures,  then 
cool  to  less  than  104  K  in  minutes  while  remaining  highly  ionized.  To  reliably  analyze  such 
data  requires  a  solid  theoretical  understanding  of  the  annihilation  process. 

Positron  production  in  solar  flares  has  been  intensively  investigated  (Ramaty  et  al.  1975; 
Kozlovsky  et  al.  1987,  2004)  and  the  positrons  were  found  to  result  from  nuclear  interactions 
of  flare-accelerated  particles  with  the  solar  atmosphere.  The  interactions  produce  radioactive 
isotopes  and,  if  there  are  significant  numbers  of  particles  with  energies  of  a  few  hundred  MeV, 
7 r+  mesons,  all  of  which  decay  to  produce  positrons.  The  initial  energies  of  the  positrons 
depend  on  the  specific  emitter:  several  hundred  keV  for  radioactive  isotopes  (e.g.,  positrons 
from  nC  and  150  have  /3-decay  spectra  with  endpoint  energies  of  ~0.97  and  1.74  MeV, 
respectively)  and  from  tens  to  hundreds  of  MeV  for  7r+  (Murphy  et  al.  1987).  In  all  cases,  the 
initial  positron  energies  are  much  higher  than  the  thermal  energies  of  the  ambient  particles. 
The  positrons  lose  energy  continuously  via  Coulomb  collisions  with  any  free  electrons  and 
discretely  via  excitation  and  ionization  of  any  neutral  atoms.  (Radiation  losses  only  become 
important  at  energies  approaching  100  MeV.) 

The  positrons  eventually  annihilate  with  electrons,  either  directly  or  through  the  forma¬ 
tion  of  a  bound  positronium  atom  (hereafter  Ps).  Direct  annihilation  of  a  positron-electron 
pair  at  rest  results  in  the  creation  of  two  y-ray  photons,  each  of  energy  e  =  mec 2  =  511  keV. 
If  Ps  is  formed,  the  annihilation  occurs  from  the  ground  state  with  the  radiative  outcome 
depending  upon  the  total  spin  wave  function  of  the  Ps  atom.  There  are  four  possible  spin 
configurations  of  the  ground  state,  one  corresponding  to  total  spin  0  (singlet)  and  three  cor¬ 
responding  to  total  spin  1  (triplet)  in  units  of  h.  Since  each  occurs  with  equal  probability,  1/4 
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of  the  Ps  atoms  will  be  observed  to  be  in  the  spin  singlet  state  (parapositronium,  hereafter 
3Ps)  and  3/4  of  the  Ps  atoms  will  be  observed  to  be  in  the  spin  triplet  state  (orthopositron- 
ium,  hereafter  3Ps).  The  singlet  state  has  a  lifetime  77  =  1.25  x  1CT10  s  and  its  annihilation 
results  in  the  creation  of  two  511  keV  photons  (27)  and  the  formation  of  a  spectral  line.  The 
triplet  state  has  a  longer  lifetime  73  =  1.42  x  10“'  s  and  its  annihilation  yields  three  photons 
(3y)  with  energies  e  <  511  keV  and  produces  a  continuum  spectrum. 

The  only  comprehensive  theoretical  study  of  positron  annihilation  in  solar  flares  was 
done  over  30  years  ago  in  the  pioneering  work  of  Crannell  et  al.  (1976).  They  identified 
for  the  first  time  the  various  annihilation  processes  relevant  for  solar  flares.  They  calculated 
the  2y  511-ken  line  width  and  the  flux  ratio  of  the  37  continuum  to  the  27  line  (hereafter, 
Q37/Q27)  in  H  plasmas,  assumed  to  be  either  neutral,  completely  ionized,  or  partially  ionized 
(10%  and  50%).  Their  calculations  included  depletion  of  3Ps  (quenching)  due  to  collisional 
interactions  with  neutral  H  (see  discussion  below). 

Subsequent  to  this  solar  work,  attention  shifted  to  the  transport  and  annihilation  of 
positrons  in  the  interstellar  medium  (ISM)  following  the  balloon  and  satellite  observations  of 
a  galactic  annihilation  line  (see  review  by  Lingenfelter  &  Ramaty  (1989)).  In  a  seminal  work, 
Bussard  et  al.  (1979)  (hereafter  BRD)  identified  the  processes  and  estimated  the  parameters 
relevant  to  the  ISM.  They  also  introduced  the  Monte-Carlo  technique  for  calculating  positron 
transport  and  Ps  formation  which  has  been  used  in  all  subsequent  work.  This  interest  in 
the  Galactic  annihilation  line  led  to  further  studies  of  annihilation  in  tenuous  media  (e.g. 
Brown  &  Leventhal  (1987);  Zurek  et  al.  (1985);  Guessoum  et  al.  (1991);  Wallyn  et  al. 
(1994);  Chapuis  et  al.  (1994);  Guessoum  et  al.  (2004)).  Recently,  Guessoum  et  al.  (2005) 
presented  updated  calculations  of  positron  annihilation  in  the  ISM,  including  the  effects  of 
molecular  H  and  grains.  Churazov  et  al.  (2005)  also  performed  new  calculations  of  positron 
annihilation  in  a  low-density  medium  and  applied  them  to  SPI/INTEGRAL  observations  of 
the  Galactic  Center  (Jean  et  al.  2003).  Such  calculations  can  be  applied  to  solar  flares, 
but  the  higher  densities  and  temperatures  often  found  in  the  solar  environment  require  the 
inclusion  of  3Ps  quenching. 

The  complex  network  of  processes  associated  with  positron  annihilation  in  solar  flares 
is  illustrated  in  Figure  1  and  the  various  interactions  are  listed  in  Table  1.  In  flares,  the 
four  relevant  annihilation  processes  are:  (1)  direct  annihilation  with  free  electrons  ( daf ),  (2) 
direct  annihilation  with  bound  electrons  (dab),  (3)  Ps  formation  via  charge  exchange  with 
neutral  atoms  (ce),  and  (4)  Ps  formation  via  radiative  combination  with  free  electrons  (rc). 
As  the  positrons  slow  down  from  their  initial  energies,  a  fraction  of  them  will  form  Ps  by 
charge  exchange  with  neutral  atoms  (process  3)  before  they  have  thermalized.  (We  note  that 
while  all  four  of  the  annihilation  processes  are  possible  before  thermalization,  only  charge 
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exchange  is  significant,  as  will  be  shown  in  §3.)  When  Ps  is  formed  in  this  way,  the  positron 
velocity  can  be  much  greater  than  typical  thermal  velocities  and  the  resulting  annihilation 
radiation  can  be  significantly  Doppler  broadened.  The  remaining  positrons  will  therm alize 
and  then  annihilate  via  one  of  the  four  processes.  In  this  case,  the  velocity  distribution  of 
the  positron  is  Maxwellian,  characterized  by  the  temperature  of  the  environment.  Unless 
the  temperature  is  quite  high  (>  105  K),  this  annihilation  radiation  will  be  less  Doppler 
broadened  than  that  from  annihilation  via  charge  exchange  before  thermalization.  The 
positron  lifetime  is  therefore  divided  into  two  epochs:  before  thermalization  (referred  to  as 
“in  flight”)  and  after  thermalization. 

Once  Ps  is  formed,  either  by  in-flight  charge  exchange  (denoted  ‘i/ce’)  or  after  thermal¬ 
ization,  depletion  of  3Ps  before  annihilation  via  collisional  quenching  can  be  significant  when 
the  density  is  sufficiently  high,  as  in  the  deep  layers  of  the  solar  atmosphere  where  positron 
production  and  subsequent  annihilation  can  take  place.  (Because  the  lifetime  of  3Ps  is  a 
thousand  times  shorter  than  that  of  3Ps,  the  required  densities  are  a  thousand  times  greater 
and  so  quenching  is  typically  only  relevant  for  3Ps.)  Quenching  is  caused  by  either  ioniza¬ 
tion  of  3Ps  or  by  conversion  of  3Ps  to  4Ps  (spin  flip).  The  density  above  which  quenching 
becomes  significant  can  be  estimated.  The  cross  section  for  the  interaction  is  approximated 
by  a  =  7TOq  (where  a q  is  the  Bohr  radius  of  the  H  atom)  and  the  quenching  rate  is  then 
rincrv,  where  Tin  is  the  total  hydrogen  number  density  and  v  is  the  relative  velocity,  approx¬ 
imated  by  a JkT /me.  For  3Ps  quenching  to  be  important,  this  rate  must  be  larger  than  the 
rate  of  annihilation  from  the  triplet  state  (l/r3),  so  that  «h  >  ( nalr^y/kT /me)~l .  At  low 
densities  and  temperatures,  such  as  those  found  in  the  ISM,  such  quenching  processes  can  be 
neglected  and  the  fraction  of  positrons  which  form  Ps  (fps)  and  Q^/Q-2^  are  directly  related 
by  Q37/Q27  =  2.25  fPs  /  (2  —  1.5  fPs). 

Quenching  cannot  be  neglected  in  the  solar  case.  In  a  plasma  of  temperature  104 
K,  relevant  to  the  solar  atmosphere,  the  above  estimate  shows  that  quenching  becomes 
important  at  densities  np  —  1015  cm-3  which  is  typical  of  the  chromosphere,  and  the  direct 
relationship  of  fps  and  1 0,2^  is  broken.  We  will  see  that  many  of  the  quenching  cross 
sections  are,  in  fact,  larger  than  7TOq  so  the  minimum  density  needed  for  quenching  can 
be  considerably  less  than  this  estimate.  In  addition,  Ps  formed  in  flight  can  have  much 
higher  velocities  relative  to  the  ambient  medium  which  also  reduces  the  density  needed  for 
quenching.  Spin-flip  quenching  reduces  the  number  of  3y  continuum  photons  and  increases 
the  number  of  2y  line  photons.  Ionization  quenching  ejects  positrons  which  can  then  form 
Ps  again  in  flight  before  they  thermalize. 

From  the  above  discussion,  positron  annihilation  in  solar  flares  is  seen  to  be  complex  and 
the  annihilation  radiation  is  correspondingly  rich  in  information  on  the  temperature,  density, 
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and  ionization  state  of  the  ambient  medium  where  the  positrons  annihilate.  Extracting 
the  most  useful  information  from  annihilation  radiation  requires  high  spectral  resolution 
measurements  such  as  those  obtained  with  RHESSI  which,  in  turn,  requires  that  the  theory 
of  annihilation  in  flares  be  updated.  We  report  here  on  the  results  of  such  a  effort  with 
a  focus  on  the  understanding  of  the  annihilation  process  after  the  positrons  slow  down  to 
energies  close  to  the  thermal  energy  in  the  local  plasma.  Detailed  treatment  of  the  transport 
and  energy  loss  of  the  high-energy  positrons  after  they  are  produced  will  be  the  subject  of  a 
separate  investigation. 

In  this  work  we  follow  the  theoretical  framework  for  positron  annihilation  developed  by 
Crannell  et  al.  (1976)  and  BRD  but  with  significant  improvements.  Since  those  pioneering 
papers  there  has  been  an  explosion  of  both  laboratory  measurements  and  calculations  of  cross 
sections  for  most  of  the  processes  involving  positrons  (see  Charlton  &  Humberston  (2001)). 
We  therefore  update  all  cross-sections  relevant  for  positron  annihilation.  We  also  extend  the 
relevant  reactions  that  were  developed  for  H  to  include  those  for  He.  Although  He/H  in  the 
Sun  is  only  ~10%,  there  can  be  regions  in  the  solar  atmosphere  where  He  remains  neutral 
but  H  is  fully  ionized  so  that  Ps  formation  via  charge  exchange  with  He  can  be  an  important 
process. 

The  calculations  have  been  performed  using  computer  algorithms  that  contain  all  the 
relevant  physics  associated  with  positron  annihilation  in  solar  flares.  We  Erst  use  an  algo¬ 
rithm  employing  a  Monte-Carlo  technique  to  calculate  (1)  the  fraction,  /l5  of  positrons  that 
form  Ps  by  in-flight  charge  exchange,  (2)  the  resulting  Ps  energy  distribution,  and  (3)  the 
3Ps  2y  line  shape.  The  remaining  positrons  thermalize,  and  we  use  several  algorithms  to 
calculate  their  interaction  rates  (da/,  dab ,  ce  and  re)  and  their  1Ps  2  7  line  shapes  as  functions 
of  the  ambient-medium  temperature.  We  then  use  several  algorithms  to  calculate  the  3Ps 
quenching  rates  as  functions  of  the  ambient-medium  temperature  for  the  3Ps  energy  distri¬ 
butions  from  each  of  the  four  Ps  formation  processes:  rc,  ce  with  H,  ce  with  He,  and  ifee. 
For  those  positrons  ejected  by  ionization  quenching  of  3Ps  formed  in-flight,  we  reapply  the 
Monte  Carlo  technique  to  calculate  (1)  the  fraction,  /2,  of  these  “second-generation”  (2 g) 
positrons  that  form  Ps  again  before  thermalization,  (2)  the  resulting  Ps  energy  distribution, 
(3)  and  the  1Ps  2 7  line  shape.  We  also  calculate  the  quenching  rate  of  this  second-generation 
Ps  but  we  assume  that  any  “third-generation”  positrons  ejected  from  ionization  quenching 
will  thermalize  since  most  of  them  will  have  energies  below  the  charge-exchange  threshold 
of  6.8  eV.  A  final  algorithm  combines  the  results  of  all  of  the  rate  algorithms  to  calculate 
the  relative  contributions  of  each  process  and  the  final  combined  shape  of  the  annihilation 
radiation  spectrum.  The  algorithms  are  made  available  as  a  separate  appendix  to  this  paper. 

In  all  of  these  calculations,  we  assume  that  the  angular  distribution  of  the  positrons  when 
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they  annihilate  or  form  Ps  is  isotropic,  which  would  certainly  be  true  once  the  positrons  have 
thermalized.  Emission  of  positrons  from  radioactive  emitters  is  isotropic  in  the  emitter  rest 
frame  and  relativistic  effects  will  not  change  this  distribution  in  the  lab  frame  even  if  the 
accelerated  particles  were  beamed  because  the  recoil  velocity  of  these  nuclei  after  production 
is  low.  Also,  in  many  cases  the  emitter  lifetimes  are  long  which  will  allow  them  time  to 
slow  down  and  isotropize.  Positrons  produced  from  pion  decay,  however,  could  be  beamed  if 
the  accelerated  particles  were  beamed.  But  since  the  only  significant  annihilation-related  in¬ 
flight  interactions  occur  at  relatively  low  energies  after  many  interactions  with  free  electrons 
or  neutral  atoms,  the  positron  angular  distribution  should  again  be  essentially  isotropic.  We 
ignore  any  possible  focussing  effects  of  transport  of  the  positrons  in  the  magnetic  loop. 

In  §2  we  describe  in  detail  the  most  recent  laboratory  measurements  and  calculations  of 
the  various  cross  sections  associated  with  positron  annihilation:  direct  annihilation,  Ps  for¬ 
mation  and  quenching,  and  ionization  and  excitation  of  the  ambient  medium  by  positrons.  In 
§3  we  describe  the  Monte  Carlo  scheme  we  use  to  treat  the  in-flight  interactions  of  positrons 
prior  to  thermalization.  In  §4  we  calculate  thermally-averaged  rates  as  a  function  of  tem¬ 
perature  for  each  of  the  annihilation  modes  of  thermalized  positrons  and  for  each  of  Ps 
quenching  interactions  using  the  updated  cross  sections.  In  §5  we  describe  the  shape  of  the 
2 7  annihilation  line  resulting  from  each  of  the  annihilation  modes  and  the  3q  continuum.  In 
§6  we  use  the  rates  calculated  in  §4  to  solve  a  system  of  continuity  equations  for  positron  and 
Ps  equilibrium  which  will  provide  the  relative  rates  for  each  of  the  annihilation  modes.  In  §7 
we  apply  the  results  to  four  specific  environments  (fully  ionized,  completely  neutral,  partially 
ionized,  and  the  solar  atmosphere  model  C  of  Vernazza  et  al.  (1981)),  calculating  the  shape 
of  the  2q  line  and  the  Q^1/Q2')  ratio  for  each.  In  §8  we  apply  the  results  to  RHESSI  data 
and  in  59  we  summarize. 


2.  Cross  Sections  for  Interactions  of  Positrons  and  Positronium 

In  this  section  we  present  the  most  recent  measurements  and  calculations  of  cross  sec¬ 
tions  for  the  interactions  relevant  to  positron  annihilation  in  solar  flares.  Uncertainties  of 
measured  cross  sections  are  typically  less  than  15%.  Uncertainties  of  calculated  cross  sec¬ 
tions  are  estimated  at  less  than  10%  (Walters,  priv.  comm.).  Where  extrapolation  to  higher 
energies  was  required,  the  energy  dependence  is  modeled  on  that  of  similar  reactions  of  elec¬ 
trons  (for  which  measurements  generally  exist)  since  the  behavior  of  electrons  and  positrons 
at  higher  energies  should  be  similar.  Cross  sections  for  direct  annihilation  of  positrons  with 
electrons  are  presented  in  §2.1,  for  Ps  formation  in  §2.2,  and  for  Ps  quenching  in  §2.3. 
Positron  energy-loss  processes  are  discussed  in  §2.4. 
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2.1.  Direct  Annihilation 


In  the  plane-wave  approximation,  the  direct  annihilation  of  a  positron  and  an  electron 
with  the  emission  of  two  photons  (e+e_  — >  27)  has  a  total  cross  section  in  the  center  of  mass 
given  by  (e.g.  Heitlcr  (1954)) 


a 


* 


^0 

4y*2/3* 


2)  + 


3  —  /3*4  ,  1  +  /3* 

- In - 

p*  1-/5* 


(i) 


Here  (3*c  is  the  magnitude  of  the  velocity  of  either  the  electron  or  positron  in  the  center  of 
mass  system,  7*  =  (1  —  /3*2)-1//2,  and  r 0  =  e2/mec2  =  2.82  x  10-13  cm.  (Hereafter,  asterisks 
refer  to  the  center  of  mass  system.) 

For  the  nonrclativistic  region  where  most  of  the  annihilations  take  place,  this  equa¬ 
tion  reduces  to  cr*  =  nr^c/v.  At  these  nonrelativistic  energies,  the  Coulomb  attraction 
between  the  particles  modifies  the  ingoing  wave  functions,  requiring  the  inclusion  of  the  fac¬ 
tor  (2na  /  (3*)  /  (l  —  e~2na^*)  (see  Crannell  et  al.  (1976)),  where  the  fine  structure  constant 
a  =  e2 /he  ~  1/137.  This  cross  section  is  shown  in  Figure  2  (multiplied  by  105  for  plotting 
purposes).  The  cross  section  is  much  less  than  the  cross  sections  for  charge  exchange  (see 
§2.2),  and  rises  with  decreasing  positron  energy. 

We  do  not  explicitly  show  the  cross  section  for  direct  annihilation  with  bound  electrons. 
The  calculation  of  the  thermally-averaged  rate  for  this  process  will  be  discussed  in  §4. 


2.2.  Positronium  Formation  Processes 

Positrons  can  form  Ps  by  charge  exchange  with  atomic  H  or  He  (e+  H  — >  H+  Ps  or  e+ 
He  — >  He+  Ps)  and  by  radiative  combination  with  free  electrons  (e+  e~  — >  7  Ps).  Positronium 
formation  through  charge  exchange  with  H  has  a  threshold  energy  of  6.8  eV.  In  Figure  2  we 
show  the  cross  section  measured  by  Zhou  et  al.  (1997)  from  threshold  to  greater  than  100 
eV  where  it  becomes  negligible.  Near  the  peak,  this  cross  section  is  significantly  less  than 
that  measured  by  Sperber  et  al.  (1992)  (used  by  Wallyn  et  al.  (1994)  and  Guessoum  et  al. 
(1997)),  but  the  Sperber  measurement  was  later  superseded  by  a  new  measurement  by  the 
same  group  (Weber  et  al.  1994)  which  is  consistent  with  the  Zhou  et  al.  (1997)  result. 

The  cross  section  for  Ps  formation  through  charge  exchange  with  neutral  He  was  mea¬ 
sured  from  threshold  (17.8  eV)  to  250  eV  by  Overton  et  al.  (1993)  and  is  also  shown  in 
Figure  2.  Above  100  eV,  this  cross  section  is  less  than  that  used  by  Chapuis  et  al.  (1994) 
who  used  earlier  measurements. 


We  do  not  explicitly  show  the  cross  section  for  Ps  formation  via  radiative  combination 
with  free  electrons.  The  calculation  of  the  rate  for  this  process  will  be  discussed  in  §4. 

We  neglect  in-flight  charge  exchange  with  He+  because  (1)  the  cross  section  for  this 
reaction  (involving  two  positively-charged  particles)  is  expected  to  be  much  smaller  than 
that  for  neutral  He,  and  (2)  in  an  environment  where  there  is  significant  He+,  H  would  be  be 
completely  ionized  and  the  increased  positron  energy  loss  due  to  free  electrons  would  favor 
thermalization  over  charge  exchange  in  flight  with  He+  (and  He/H  is  only  0.1).  We  also 
neglect  charge  exchange  of  thermalized  positrons  with  He+  because  (1)  of  its  significantly 
higher  threshold  energy  (47.6  eV)  and  (2)  the  ionized  environment  would  favor  Ps  formation 
by  radiative  combination  with  free  electrons. 


2.3.  Positronium  Quenching  Processes 

The  cross  sections  for  ionization  and  spin  flip  of  3Ps  by  free  electrons  ( e~  3Ps  — >  e~e~e+ 
and  e-  3Ps  — >  e~  3Ps)  were  recently  calculated  up  to  80  eV  by  Gilmore  et  al.  (2004)  and 
are  shown  in  Figure  3.  This  is  an  important  interaction  in  ionized  gases  at  densities  >  1013 
cm-3,  typical  of  the  chromosphere.  In  their  treatment  of  3Ps  quenching,  Crannell  et  al. 
(1976)  included  only  interactions  with  neutral  H  and  hence  underestimated  the  effects  of 
quenching. 

Ionization  of  3Ps  by  neutral  H  (H  3Ps  — >  H  e_e+)  was  included  in  the  analysis  of  Crannell 
et  al.  (1976)  where  the  cross  section  calculated  in  the  Born  approximation  by  Massey  & 
Mohr  (1954)  was  used.  Here,  we  use  more  sophisticated  calculations,  using  coupled  pseudo¬ 
state  formalism,  recently  performed  by  Campbell  et  al.  (1998b)  up  to  40  eV.  Above  that 
energy  we  use  the  calculations  of  Ray  et  al.  (2002).  The  resulting  cross  section  is  shown  in 
Figure  3.  The  cross  section  for  spin  flip  of  3Ps  by  H  (H  3Ps  — >  H  3Ps)  was  recently  calculated 
by  Blackwood  et  al.  (2002)  and  is  also  shown  in  Figure  3. 

Much  effort  has  gone  into  obtaining  the  cross  section  for  the  inverse  process  of  charge 
exchange  with  ionized  H  (H+  Ps  — >  H  e+)  since  the  charge  conjugate  reaction  p  Ps  — > 
He-  results  in  the  formation  of  anti-hydrogen.  Because  of  time  reversal  invariance,  the 
cross  section  for  inverse  charge  exchange  is  related  to  that  for  charge  exchange  through  (e.g. 
Humberston  et  al.  (1987)) 

E 

u(H+  Ps  ->  H  e+)  =  — ^n(H  e+  -»•  H+  Ps)  (2) 

2Eps 

where  E+  and  Eps  are  the  kinetic  energies  of  the  positron  and  Ps,  respectively.  Humberston 
et  al.  (1987)  derived  a  cross  section  for  this  process  from  theoretical  calculations  of  the 
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charge  exchange  cross  section.  Here  we  use  the  measurements  of  the  charge  exchange  cross 
section  described  above  to  obtain  <v(H+  Ps  — >  H  e+).  For  inverse  charge  exchange  in  singly- 
ionized  He  (He+  Ps  — >  He  e+),  we  proceed  as  in  the  case  for  H  and  invert  the  He  charge 
exchange  cross  section  shown  in  Figure  2. 

We  neglect  inverse  charge  exchange  interactions  with  doubly-ionized  He.  The  foremost 
reason  for  this  is  that  if  there  is  any  appreciable  ionization  then  free  electrons  will  dominate 
the  quenching.  Also,  the  presence  of  doubly-ionized  He  would  imply  high  temperature  where 
direct  annihilation  (with  free  electrons)  would  dominate  over  Ps  formation  (by  radiative 
combination).  This,  coupled  with  the  fact  that  He  is  lower  in  abundance  than  H  by  an  order 
of  magnitude,  justifies  neglecting  quenching  by  doubly-ionized  He. 

We  neglect  spin-flip  and  ionization  quenching  of  Ps  by  neutral  He  because  (1)  if  the 
medium  is  neutral,  ionization  by  neutral  H  will  dominate  since  He/H  =  0.1  and  (2)  if  the 
medium  is  ionized,  quenching  by  free  electrons  will  dominate.  In  addition,  spin-flip  interac¬ 
tions  resulting  from  electron  exchange  between  Ps  and  neutral  He  have  a  high  threshold  of 
19.8  eV.  This  high  threshold  results  from  the  fact  that  in  the  ground  state  the  electrons  in 
the  He  atom  must  form  a  spin  singlet  and  electron  exchange  interactions  which  convert  3Ps 
to  3Ps  require  exciting  He  to  the  triplet  state  (Fraser  &  Krady  1965).  However,  “pick-off” 
quenching  can  occur,  where  the  Ps  positron  annihilates  directly  with  a  bound  electron  in 
which  it  shares  a  relative  spin  singlet  state.  The  experimentally  determined  rate  is  pressure 
dependent  (Fraser  &  Krady  1965).  Using  the  experimental  results  of  Beers  &  Hughs  (1968) 
and  the  ideal  gas  law,  the  quenching  rate  is  7.3  x  10~20nT  s^1  with  n  in  cm-3.  For  this  to 
exceed  the  rate  for  annihilation  requires  nT  >  1026  K  cm-3  which  is  unlikely  in  the  flare 
plasma. 

We  also  neglect  collisional  spin  flip  due  to  interactions  with  ions  because  the  cross 
sections  are  expected  to  be  significantly  lower  than  those  for  interactions  where  electrons  are 
exchanged.  Also,  we  do  not  include  spin  flip  of  3Ps  due  to  magnetic  fields.  However,  as  has 
been  measured  by  Pond  &  Dicke  (1952),  there  can  be  additional  quenching  by  as  much  as 
30%  for  magnetic  held  strengths  of  several  thousand  gauss. 


2.4.  Positron  Energy  Loss 

Positrons  can  lose  energy  via  discrete  interactions  with  the  neutral  component  of  the 
plasma.  The  energy  lost  by  the  positron  in  these  interactions  is  discussed  in  §3.1.  Here 
we  discuss  the  interaction  cross  sections.  In  Figure  2  we  present  the  cross  sections  for  the 
ionization  reactions  e+  H  — >  e+e_  H+  and  e+  He  — >  e+e”  He+.  (These  reactions  are  typically 
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referred  to  as  “impact”  ionization  to  distinguish  them  from  ionization  associated  with  Ps 
formation  through  charge  exchange,  e+  H  — >  H+  Ps  and  e+  He  — »  He+  Ps,  which  were 
treated  in  §2.2.)  The  cross  section  for  impact  ionization  of  H  was  measured  by  Jones  et  al. 
(1993)  from  threshold  to  700  eV  and  confirmed  by  Hoffman  et  al.  (1997).  These  results 
are  significantly  lower  than  those  of  Spicher  et  al.  (1990)  which  were  used  by  Wallyn  et  al. 
(1994)  in  their  study  of  annihilation  in  the  ISM.  The  cross  section  for  impact  ionization  of 
He  was  measured  by  Jacobsen  et  al.  (1995)  from  threshold  to  1500  eV.  From  threshold  to 
~100  eV,  these  results  are  somewhat  less  than  earlier  published  data  (e.g.,  those  of  Fromme 
et  al.  (1986)  used  by  Chapuis  et  al.  (1994)  who  first  included  He  reactions  for  the  ISM). 

There  are  no  measurements  of  positron  excitation  of  H  (i.e.,  e+  H  — »  e+  H*).  However, 
the  elaborate  calculations  of  Kernoghan  et  al.  (1996)  successfully  reproduced  the  available 
data  for  the  total,  elastic,  ionization  and  positronium-formation  cross  sections.  This  success 
lends  assurance  that  their  calculation  of  the  ls-2s  and  ls-2p  excitation  cross  sections  are 
reliable  as  well.  Stein  et  al.  (1998)  estimated  the  cross  section  for  exciting  the  higher 
states  of  H  by  subtracting  the  partial  cross  sections  for  clastic,  ionization,  ls-2s  and  ls-2p 
excitation,  and  Ps  formation  from  the  total  cross  section.  In  Figure  2  we  show  the  resulting 
cross  section  for  excitation  to  all  levels  of  H. 

Mori  et  al.  (1994)  measured  the  total  cross  section  for  excitation  of  He  by  positrons 
(e+  He  — >  e+  He*)  and  it  is  shown  in  Figure  2.  The  recent  calculation  of  Wu  et  al.  (2004) 
is  in  general  agreement  with  these  data.  We  note  that  these  results  are  significantly  less 
than  the  cross  section  evaluated  by  Chapuis  et  al.  (1994);  their  evaluation  assumed  that 
elastic  scattering  is  negligible,  as  it  is  for  H.  But,  as  discussed  by  Campbell  et  al.  (1998a), 
clastic  scattering  is  significant  for  He  and  ionization  is  the  largest  cross  section  above  75  eV, 
in  contrast  to  the  result  of  Chapuis  et  al.  (1994)  in  which  excitation  had  the  largest  cross 
section. 

We  do  not  include  excitation  and  ionization  of  He+  because  if  He+  is  abundant,  H  will 
be  mostly  ionized  and  positron  energy  loss  will  be  dominated  by  Coulomb  interactions  with 
free  electrons. 


Positrons  can  also  lose  energy  continuously  to  the  ionized  component  of  the  plasma. 
The  energy-loss  rate  is  given  by  the  expression  (Huba  1994) 


dE 

dt 


1.7  x  10 ~*ne(3~l  In  A 


( E/kT)1/2e-E,iT 


(eV/s),  (3) 


where  E  and  /3c  are  the  positron  energy  (eV)  and  velocity  and  T  and  ne  are  the  temperature 
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3.  In-Flight  Positronium  Formation 

When  first  produced,  positrons  have  initial  energies  of  ~1  MeV  (from  decay  of  radioac¬ 
tive  positron  emitters)  to  hundreds  of  MeV  (from  decay  of  pions).  These  positrons  lose 
energy  to  the  ambient  medium  via  Coulomb  interactions  with  free  electrons  and  via  ioniza¬ 
tion  and  excitation  of  neutral  atoms.  Coulomb  interactions  with  free  electrons  will  always 
dominate  over  direct  annihilation  or  radiative  combination  with  free  electrons  because  these 
latter  processes  involve  the  emission  of  photons.  The  cross  section  for  direct  annihilation 
with  bound  electrons  at  high  energies  is  simply  the  cross  section  for  annihilation  with  free 
electrons  (equation  1)  which,  at  a  few  hundred  keV,  is  on  the  order  of  10~26  cm2.  Extrapolat¬ 
ing  the  cross  section  for  H  ionization  by  positrons  of  Figure  2  to  these  energies  as  E-1  (Allen 
1976)  gives  ICC21  cm2.  Positrons  at  these  energies  will  therefore  lose  energy  via  ionization 
rather  than  annihilating  with  free  or  bound  electrons.  The  small  fraction  of  positrons  that 
do  annihilate  at  higher  energies  will  produce  an  annihilation  line  that  will  be  significantly 
Doppler  broadened  and  buried  in  the  underlying  continuum.  The  only  remaining  process  is 
charge  exchange  whose  cross  section  becomes  significant  at  energies  less  than  ~100  eV  and 
dominates  the  other  processes. 

Positrons  can  therefore  form  Ps  through  charge  exchange  with  neutral  H  and  He  as  they 
slow  down  if  the  ambient  medium  is  not  fully  ionized.  At  densities  greater  than  ~1013  cm-3, 
quenching  via  Ps  ionization  results  in  ejected  positrons  that  can  again  form  Ps  in  flight. 
We  refer  to  these  latter  positrons  as  “second-generation”  (2 g)  to  distinguish  them  from  the 
“primary”  positrons.  In  §3.1  we  discuss  our  calculations  of  in-flight  Ps  formation  by  primary 
positrons  and  in  §3.2  we  discuss  in-flight  Ps  formation  by  second-generation  positrons. 

3.1.  Positronium  Formation  by  Primary  Positrons 

We  use  a  Monte  Carlo  simulation  to  calculate  the  fraction  of  positrons  that  form  Ps 
prior  to  thermalization  via  charge  exchange  (/i),  the  energy  distribution  of  the  resulting  Ps, 


and  free  electron  density  (cm  3).  The  Coulomb  logarithm  is  given  by 
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and  the  shape  of  the  27  line  from  the  decay  of  1Ps  in  flight.  The  procedure  is  similar  to  that 
employed  by  BRD  except  that  we  include  He  and  use  the  updated  cross  sections  discussed 
above.  In  the  simulation,  positrons  are  injected  into  the  plasma  with  an  initial  energy  of 
700  eV.  The  results  are  not  sensitive  to  the  exact  initial  energy  as  long  as  it  is  sufficiently 
above  energies  in  the  range  10-100  eV  where  charge  exchange  is  important.  The  positrons 
lose  energy  continuously  to  the  ionized  component  of  the  plasma  until  they  interact  with  a 
neutral  H  or  He  atom  where  they  can  excite,  ionize  or  capture  a  bound  electron  to  form  Ps. 


Given  an  initial  energy  Eq  and  a  uniformly-distributed  random  number  r  from  the  unit 
interval,  the  energy  E  at  which  the  positron  interacts  discretely  is  determined  from 


r  =  exp 


— cnH 


dE 


fob 


Ot 


dE/dt\ 


(6) 


where  dE/dt  is  given  by  equation  (3)  and  atot  is  the  total  cross  section  for  making  any  kind 
of  discrete  interaction  with  target  H  or  He  atoms,  given  by 


0tot  =  (1  -  XH  +)(<t®  +  afon  +  ofj  +  (nHe/nH)  (1  -  ^He+  -  XRc ++)(a*e  +  +  acHee).  (7) 


where  X,  is  the  ratio  of  the  number  of  ions  of  species  i  to  the  total  number  of  atoms  of  the 
given  element  (ionized  plus  neutral).  The  probability  that  a  given  type  of  discrete  interaction 
occurs  is  then  proportional  to  the  cross  section  for  that  interaction.  The  cross  sections  were 
discussed  in  §2  and  are  shown  in  Figure  2. 

The  energy  lost  by  the  positron  in  interactions  with  neutral  atoms  is  calculated  as 
follows.  As  in  BRD,  in  collisions  which  excite  H  the  average  amount  of  energy  lost  by  the 
positron  is  taken  to  be  10.2  eV,  the  excitation  energy  of  the  2xs  level.  For  ionization  collisions 
with  H  the  ejected  electron  velocity  distribution  is  approximately  Gaussian  with  a  standard 
deviation  of  0.382ac  ~  2  eV  (where  a  is  the  fine  structure  constant)  and  centered  at  one 
fourth  of  the  binding  energy  (BRD;  Omidvar  (1965)).  Conservation  of  energy  then  provides 
the  positron  energy.  For  He  we  follow  Chapuis  et  al.  (1994)  and  adopt  a  mean  energy  loss 
per  excitation  of  20.6  eV,  which  corresponds  to  excitation  from  the  ground  state  to  the  2xs 
level,  and  we  use  the  value  of  49  eV  (Opal  et  al.  1971)  for  the  mean  energy  lost  in  an 
ionization  collision. 

After  reducing  the  positron  energy  according  to  these  prescriptions,  the  algorithm  rein¬ 
jects  the  positron  and  the  above  process  is  repeated  until  one  of  the  three  following  outcomes 
occurs:  (1)  the  energy  of  the  positron  falls  below  6.8  eV,  the  threshold  for  charge  exchange 
with  H;  (2)  the  positron  reaches  thermal  energies  defined  by  dE/dt  >  0,  i.e.  the  positron  no 
longer  loses  energy;  or  (3)  a  charge-exchange  interaction  occurs.  If  either  (1)  or  (2)  happens 
then  the  positron  is  recorded  as  thermalizing.  If  (3)  occurs  then  the  charge  exchange  event 
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is  counted  and  the  Ps  energy  is  calculated,  recorded  and  binned,  where  the  energy  of  the  Ps 
is  the  current  energy  of  the  positron  less  the  difference  between  the  ionization  energy  of  H 
(13.6  eV)  or  He  (24.6  eV  )  and  the  Ps  binding  energy  (6.8  eV).  The  entire  process  is  then 
repeated  until  sufficiently  good  statistics  are  obtained  for  the  Ps  energy  distribution. 

In  Figure  4  we  show  the  fraction  f\  of  positrons  that  form  Ps  in  flight  (solid  curve)  as 
a  function  of  the  ionization  fraction  of  H,  A"H+.  The  temperature  and  electron  density  are 
fixed  at  T=ll,600  K  and  ne  —  5  x  1013  cm-3  and  He/H=0,  values  assumed  by  BRD  in  their 
calculation  of  /j.  For  a  completely  neutral  medium  we  find  f\  =  0.95  which  is  the  same  as 
the  value  0.95  found  by  BRD.  The  effect  of  inclusion  of  He  (He/H  =  0.1)  is  shown  by  the 
dotted  curve.  The  dot-dashed  curve  is  discussed  in  §3.2. 

Our  results  for  f\  are  consistently  higher  than  those  of  BRD  (dashed  line  in  Figure  4) 
for  Xh+  >  0.05  even  though  our  charge  exchange  cross  section  is  very  similar,  our  ionization 
and  excitation  cross  sections  are  larger  (which  should  reduce  /i),  and  we  use  an  identical 
electron  energy  loss  expression.  To  confirm  our  result,  we  have  developed  a  second  Monte 
Carlo  code  that  directly  follows  positrons  with  initial  energy  of  700  eV  through  small,  finite 
path  intervals  A/.  For  each  interval  we  calculate  the  probabilities  for  each  of  the  three  discrete 
reactions  (%  &i(E+)  A l  where  i  is  ce,  ionization  or  excitation)  at  the  current  positron  energy 
( E+ ).  If  charge  exchange  occurs,  the  event  is  counted.  If  ionization  or  excitation  occurs,  the 
energy  is  reduced  according  to  the  above  prescriptions  and  the  next  interval  is  considered.  If 
no  discrete  interaction  occurs,  the  positron  energy  is  reduced  by  A E  =  ( dE/dl )  A l  (where 
dE/dl  =  ( dE/dt)v~l )  and  the  next  interval  is  considered,  unless  the  positron  energy  has 
fallen  below  the  6.8  eV  threshold  or  the  positron  has  reached  thermal  energies  defined  by 
dE/dt  >  0.  This  process  is  repeated  until  reasonably  good  statistics  are  obtained.  We  find 
that  our  new  code  reproduces  our  original  results  for  fi  exactly.  We  have  also  compared 
our  results  with  the  recent  calculations  of  positron  annihilation  in  the  interstellar  medium 
by  Guessoum  et  al.  (2005)  and  find  excellent  agreement.  We  have  no  explanation  for  the 
differences  with  BRD,  but  we  note  that  if  we  increase  our  dE/dt  by  a  factor  2  independent 
of  energy,  we  reproduce  the  BRD  results  for  f\ . 

The  solid  histogram  in  Figure  5  shows  the  energy  distribution  of  Ps  formed  in  flight  (Fps) 
in  a  neutral  medium  (He/H  =  0.1).  There  is  only  a  weak  dependence  of  this  distribution  on 
temperature  since  most  of  the  Ps  is  produced  in  the  narrow  range  of  energies  between  the 
threshold  for  charge  exchange  and  excitation  of  H  (the  Ore  gap  (Charlton  &  Humberston 
2001);  see  Figure  2)  where  the  only  significant  process  is  charge  exchange.  If  the  medium  is 
only  slightly  ionized,  Coulomb  energy  losses  rapidly  thermalize  positrons  just  above  the  Ps 
formation  threshold.  This  removes  some  of  the  lowest-energy  Ps,  slightly  modifying  the  Ps 
formation  energy  distribution  as  shown  by  the  dotted  histogram  in  Figure  5  for  Ah+  =  0.05 
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(A"He+  —  -AHe++  —  0). 

From  equation  (3),  the  estimated  energy-loss  time  (E  -y  dE/dt )  of  a  30  eV  positron  is 
>  10-6  s  at  an  electron  density  of  1012  cm-3  (corresponding  to  an  H  density  of  1014  cm-3  and 
1%  ionization  fraction).  Because  Ps  is  both  neutral  and  heavier  than  electrons,  the  energy- 
loss  time  of  30  eV  Ps  will  be  much  longer  than  this  and  that  of  higher-energy  Ps  even  longer. 
Since  these  times  are  longer  than  the  lifetimes  of  either  :Ps  or  3Ps,  Ps  formed  in-flight  will 
not  lose  significant  energy  before  annihilating  or  undergoing  quenching.  Therefore,  both  the 
Doppler-broadened  shape  of  the  2y  line  from  4Ps  decay  and  the  3Ps  quenching  rate  will  be 
determined  by  the  energy  distribution  of  the  Ps  at  formation  as  calculated  above.  Because  the 
velocities  of  Ps  formed  in-flight  are  larger  than  those  for  the  thermal  annihilation  processes 
(see  §5),  the  2y  line  from  :Ps  formed  in  flight  will  be  correspondingly  broader.  Quenching 
will  be  discussed  in  54. 


3.2.  Positronium  Formation  by  Second-Generation  Positrons 

Because  the  energy  distribution  of  the  Ps  at  formation  is  unchanged  before  undergo¬ 
ing  any  quenching,  positrons  liberated  via  Ps  ionization  will  still  have  energies  significantly 
above  thermal.  As  these  second-generation  positrons  lose  energy  they  may  again  form  Ps  via 
charge  exchange  before  thermalization.  We  use  the  Monte  Carlo  simulation  for  Ps  formation 
from  primary  positrons  described  above  to  also  calculate  the  fraction  of  second-generation 
positrons  that  form  Ps  via  charge  exchange  prior  to  thermalization  (/2),  the  energy  distri¬ 
bution  of  the  resulting  second-generation  Ps,  and  the  shape  of  the  2y  line  from  the  decay 
of  4Ps  in  flight.  In  the  simulation,  a  positronium  energy  is  randomly  selected  from  the 
distribution  of  Ps  energies  at  the  time  of  ionization  (Fps,ion),  an  ejected  positron  energy  is 
determined  from  this  Ps  energy,  and  the  positron  is  injected  with  this  initial  energy  into  the 
above  Monte  Carlo.  This  is  repeated  until  adequate  counting  statistics  are  achieved.  The 
distribution  FPs40n  is  obtained  from  the  proportionality: 

Aps^on  (X  FPs  V  <73^q,  (8) 

where  Fps  is  the  primary  Ps  energy  distribution  calculated  above  (§3.1),  v  is  the  relative 
velocity  (taken  to  be  the  velocity  of  the  Ps  since  it  is  generally  much  higher  than  the  thermal 
velocities  of  the  ambient  particles),  and  is  the  total  cross  section  for  Ps  breakup,  given 
by 

03^0  =  (!  -  Ah+V^o  +  (jie/nn)  +  XH+  crfnvce  +  (nHe/nH)  XHe+a*evce.  (9) 

Here,  u|Q0  is  the  cross  section  for  Ps  ionization  by  neutral  H,  crf^0  is  the  cross  section  for 
Ps  ionization  by  free  electrons,  and  cr^vce  and  cr?®  ce  are  the  cross  sections  for  inverse  charge 
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exchange  by  singly-ionized  H  and  He,  respectively.  These  quenching  cross  sections  were 
discussed  in  §2.3. 

We  estimate  the  energy  of  the  ejected  positron  in  the  following  way.  Because  most  of  the 
energies  of  the  primary  in-flight  Ps  are  above  10  eV  (see  Figure  5),  Ps  breakup  is  dominated 
by  ionization  via  either  neutral  H  or  free  electrons  (see  Figure  3).  For  these  reactions,  there 
are  three  particles  in  the  exit  channel  and  the  energy  distribution  of  the  ejected  positron  is 
expected  to  be  distributed  up  to  a  maximum  energy  Emax  with  a  peak  at  Emax/ 2,  where 
Emax  =  Eps  —  6.8  eV  (the  binding  energy  of  Ps)  and  Eps  is  the  energy  of  the  Ps  at  the  time 
of  ionization.  We  approximate  this  positron  energy  distribution  with  a  cosine  function.  We 
find  that  the  shape  of  the  resulting  2y  line  is  insensitive  to  the  exact  distribution  used. 

The  dot-dashed  curve  in  Figure  4  shows  the  fraction  /2  of  second-generation  positrons 
that  form  Ps  in  flight  as  a  function  of  A"h+,  the  ionization  fraction  of  H.  As  before,  the 
temperature  and  electron  density  are  fixed  at  T=ll,600  K  and  ne  —  5  x  1013  cm-3  and 
He/H  =  0.  For  a  completely  neutral  medium  we  find  /2  =  0.64.  The  dashed  histogram 
in  Figure  5  shows  the  energy  distribution  of  second-generation  Ps  formed  in  flight  in  a 
neutral  medium  (He/H  =  0).  Similar  to  primary  Ps,  there  is  only  a  weak  dependence  of 
this  distribution  on  temperature  since  most  of  the  Ps  is  produced  in  the  Ore  gap  (see  Figure 
2).  The  2q  line  from  annihilation  of  3Ps  produced  by  these  second-generation  positrons  will 
be  narrower  than  that  from  in-flight  Ps  formed  from  first-generation  positrons  because  the 
second-generation  Ps  energy  distribution  is  shifted  to  lower  energies,  as  can  be  seen  in  Figure 
5  (also  see  §5).  Positrons  from  ionization  quenching  of  second-generation  Ps  are  assumed  to 
be  thermalized  since  their  energies  will  generally  be  below  the  charge  exchange  threshold. 


4.  Interaction  Rates  and  Positronium  Energy  Distributions 

In  §6  we  will  develop  relationships  among  the  instantaneous  numbers  of  thermalized 
positrons  and  singlet  and  triplet  Ps  atoms  (formed  via  the  various  modes)  in  order  to  cal¬ 
culate  the  relative  steady-state  numbers  of  these  species.  These  relationships  will  involve 
thermally-averaged  rates  for  the  interactions  that  thermalized  positrons  and  Ps  and  in-flight 
Ps  experience  with  the  ambient  thermal  environment.  For  thermalized  positrons,  the  rates 
are  for  direct  annihilation  (with  both  free  and  bound  electrons)  and  Ps  formation  (via  ra¬ 
diative  combination  with  free  electrons  and  Ps  formation  via  charge  exchange  with  H  and 
He).  For  Ps,  the  rates  are  for  the  six  Ps  quenching  reactions  associated  with  each  of  the 
four  Ps  formation  modes  (radiative  combination  of  thermal  positrons  with  free  electrons, 
charge  exchange  of  thermal  positrons  with  H  and  He,  and  in-flight  charge  exchange).  In  this 
section  we  calculate  these  rates  and  the  energy  distribution  of  the  Ps  associated  with  each 
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the  Ps-producing  processes.  (The  energy  distribution  of  Ps  formed  in  flight  was  discussed 
in  §3.)  The  Ps  energy  distributions  will  be  used  in  §5  to  calculate  the  associated  4Ps  2y 
line  shapes.  In  the  following  calculations,  we  assume  the  ambient  plasma  is  comprised  of 
the  following:  e~  -  free  electrons;  H  -  neutral  hydrogen;  H+  -  ionized  hydrogen;  He  -  neutral 
helium;  He+  -  singly  ionized  helium;  and  He++  -  fully  ionized  helium. 


The  average  rate  per  unit  target  density  for  a  binary  interaction  in  a  plasma  is  defined 
by  A  =  (cry),  where  a  is  the  total  cross  section  for  the  given  process  and  the  relative  velocity 
v  =  |vi  —  V2I,  and  Vi  and  V2  are  the  velocities  of  the  interacting  particles.  The  averaging  is 
done  over  the  velocity  distributions  of  both  interacting  particles.  For  the  ambient  particle, 
this  distribution  is  always  Maxwellian  characterized  by  the  temperature.  For  reactions  of 
therm ali zed  positrons,  the  positron  distribution  is  also  Maxwellian  characterized  by  the  same 
temperature.  In  this  case,  the  interaction  rate  (for  kT  <C  mec 2)  is 


A 


v  27 rkT  J 


d3v  1  d3v2  exp 


(  mivj\ 
V  2kT  J 


exp 


(  rn2vl\ 

V  2kT  J 


va(v  1,  v2). 


(10) 


For  direct  annihilation  of  thermal  positrons  with  free  electrons,  equation  (10)  was  nu¬ 
merically  integrated  by  Crannell  et  al.  (1976)  and  Bussard  et  al.  (1979)  using  the  cross 
section  given  in  §2.1.  The  rate  calculated  by  Bussard  et  ah  (1979)  is  shown  in  Figure  6  as 
a  function  of  temperature. 

For  direct  annihilation  of  thermal  positrons  with  bound  electrons,  the  rate  of  reaction 
is  often  expressed  in  terms  of  the  effective  number  of  protons  ( Zeff )  multiplied  by  the  free 
annihilation  rate, 

A  =  7T  rlcZeff.  (11) 

Ze//,  which  varies  with  the  speed  of  the  positron,  is  a  measure  of  the  probability  of  the 
positron  being  at  the  same  position  as  one  of  the  target  electrons  and  is  calculated  from  the 
clastic  scattering  wave  function  for  the  positron-target  system.  For  H,  we  use  the  rate  given 
in  BRD  who  used  the  thermally-averaged  Zeff  calculated  by  Bhatia  et  ah  (1977).  The 
resulting  rate  is  shown  in  Figure  6.  The  curve  stops  at  104  K  since  H  will  be  mostly  ionized 
above  that  temperature.  For  He,  we  integrate  the  Zeff  given  by  Humberston  (1974)  over 
the  positron  Maxwellian  distribution  and  the  resulting  rate  is  also  shown  in  Figure  6. 

For  Ps  formation  via  radiative  combination  of  thermal  positrons  with  free  electrons  we 
use  the  rate  given  by  Gould  (1989)  with  the  Gaunt  factor  set  equal  to  0.9.  This  is  also 
shown  as  a  function  of  temperature  in  Figure  6  where  it  is  seen  that  this  process  is  more 
important  than  direct  annihilation  up  to  T  ~  7  x  105  K;  above  this  temperature,  direct 
annihilation  dominates.  The  energy  distribution  of  the  resulting  Ps  is  Maxwellian  (Gould 
1989).  The  reason  is  that  for  kT  <C  mec 2  the  added  velocity  kick  in  the  center  of  mass  of 
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the  e+e_  system  due  to  the  kinetic  energy  of  the  e+e-  pair  and  the  release  of  binding  energy 
can  be  shown  to  go  as  f3kiCk  —  (6.8  eV  +  kT)/2mec2,  whereas  the  velocity  of  the  center  of 
mass  system  goes  as  /3cm  ~  \] kT/mec 2,  hence  /3kiCk  -C  /3cm. 

For  Ps  formation  via  charge  exchange  of  thermal  positrons  with  ff  and  He,  we  integrated 
equation  (10)  using  the  cross  sections  given  in  §2.2  and  the  rates  are  also  shown  in  Figure 
6.  The  energy  distribution  of  the  resulting  Ps  is  not  Maxwellian  and  is  derived  as  follows. 
Given  a  positron  energy  E,  the  conditional  positronium  energy  distribution  is 

FPs(EPs\E)  oc  ace(E )  v(E)  5(EPs  —  (E  —  6.8 eV))  (12) 

where  we  have  assumed  that  the  energy  of  the  Ps  is  the  initial  energy  of  the  positron  less 
the  difference  between  the  ionization  energy  of  H  (13.6  eV)  and  the  Ps  binding  energy  (6.8 
eV).  For  He,  6.8  eV  is  replaced  by  17.8  eV  in  equation  (12)  since  the  ionization  energy  of 
He  is  24.6  eV  (24.6  —  6.8  =  17.8).  Integrating  equation  (12)  over  the  positron  Maxwellian 
velocity  distribution  gives  the  Ps  energy  distributions  for  H  and  He: 


Fps(E) 


Fps(E) 


1  dNPs 
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1  d,NPs 
NPs  dE 


(E)  oc  crfe(E+ 6.8eV)  v(E+ 6.8eV)  V E  +  6.8 eV  exp[— (E  +  6.8eV)/kT] 

(13) 

(E)  oc  <x*e(£+17.8eV)  v(E+17.8eV)  VE  +  17.8 eV  exp[-(£+17.8eV)/fcT]. 

(14) 


We  now  calculate  the  six  3Ps  quenching  rates  for  each  of  the  four  Ps  production  processes 
(radiative  combination  of  thermal  positrons  with  free  electrons,  charge  exchange  of  thermal 
positrons  with  H  and  He,  and  charge  exchange  in  flight).  The  six  quenching  cross  sections 
were  discussed  in  §2.3  and  are  shown  in  Figure  3.  For  each  of  the  four  Ps-formation  cases, 
the  velocity  distribution  of  the  ambient  particle  (electron,  H,  H+  or  He+)  is  Maxwellian  but 
the  velocity  distribution  of  the  Ps  depends  on  the  formation  process  as  discussed  above. 
For  Ps  formed  via  radiative  combination  of  thermal  positrons  with  free  electrons,  the  Ps 
distribution  is  also  Maxwellian  and  equation  (10)  may  be  used.  The  resulting  rates  are 
shown  in  panel  (a)  of  Figure  7.  For  Ps  formed  via  charge  exchange  of  thermal  positrons 
with  ambient  H  and  He,  the  Ps  velocity  distribution  is  given  by  equations  (13)  and  (14)  for 
H  and  He,  respectively.  Averaging  these  distributions  with  the  Maxwellian  ambient-particle 
distributions  results  in  the  quenching  rates  shown  in  panels  (b)  and  (c)  of  Figure  7.  For 
charge  exchange  in  flight,  the  primary  Ps  distribution  is  similar  to  that  shown  in  Figure  5 
and  the  resulting  quenching  rates  for  Ps  formed  from  primary  positrons  are  shown  in  panel 
d  of  Figure  7.  We  note  that  all  of  these  in-flight  quenching  rates  are  nearly  independent 
of  temperature.  This  is  because  the  velocity  of  the  Ps  (energies  of  tens  of  eV,  Figure  5)  is 
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much  higher  than  the  thermal  velocities  of  the  ambient  particles  so  that  the  relative  velocity 
is  essentially  the  velocity  of  the  Ps.  The  relative  velocity  (and  therefore  the  rate)  is  then 
independent  of  temperature.  Also,  the  calculations  for  the  in-flight  rates  are  only  performed 
up  to  2  x  105  K.  Above  this  temperature  there  will  probably  not  be  sufficient  numbers  of 
neutral  atoms  for  charge  exchange  to  be  important. 


5.  Annihilation  Radiation  Spectra 

The  shape  and  corresponding  width  of  the  positron  annihilation  line  is  different  for  each 
of  the  modes  of  annihilation  shown  in  Figure  1:  direct  with  free  and  bound  electrons  and  via 
Ps  formed  by  radiative  combination  and  by  charge  exchange  (either  after  thermalization  or 
in  flight).  The  shape  of  the  observed  total  annihilation  radiation  spectrum  (both  2y  line  and 
3y  continuum)  is  a  composite  of  these  components  and  their  relative  contributions  depend 
on  the  temperature,  density,  ionization  and  composition  of  the  ambient  medium  (as  will  be 
discussed  in  §6).  In  this  section  we  discuss  the  shapes  of  each  of  the  various  components: 
the  2y  line  is  discussed  in  §5.1  and  the  3y  continuum  in  §5.2. 


5.1.  The  Annihilation  Line 

The  shape  of  the  2y  line  resulting  from  the  direct  annihilation  of  nonrclativistic  thermal 
positrons  and  free  electrons  is  Gaussian  with  a  full- width  at  half  maximum  (FWHM)  given 
by 

FWHM  =  mc^ In  2^)  =  1.1  keV  (15) 

as  derived  in  Appendix  A.  Since  the  centers  of  the  positron  and  electron  energy  distributions 
in  the  center  of  mass  frame  are  on  the  order  of  kT ,  the  line  center  will  be  blue-shifted  to 
e  =  mec 2  +  (f)kT  where  4>  is  on  the  order  of  unity  (Gould  1989).  In  the  Sun,  the  annihilation 
line  could  therefore  be  blue  shifted  by  ~1  eV. 

Similarly,  the  shape  of  the  annihilation  line  resulting  from  the  2y  decay  of  xPs  formed 
via  radiative  combination  of  thermal  positrons  with  free  electrons  also  has  a  width  given 
by  equation  (15).  As  shown  above  (§4),  Ps  formed  in  this  way  has  a  Maxwellian  velocity 
distribution.  Since  the  photon  is  emitted  isotropically  in  the  rest  frame  of  the  Ps  and 
has  energy  in  the  lab  frame  given  by  equation  (A3),  the  photon  is  distributed  in  energy 
according  to  equation  (A10)  from  which  the  width  given  by  equation  (15)  is  obtained  (see 
also  an  alternative  derivation  by  Gould  (1989)). 
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However,  the  velocity  distribution  of  Ps  formed  through  thermal  charge  exchange  is  not 
Maxwellian  and  the  shape  of  the  2y  line  from  1Ps  decay  is  not  Gaussian.  To  derive  the 
line  shape  for  this  decay  (xPs  — >  2y),  we  follow  the  derivation  by  Stecker  (1971)  for  the 
analogous  decay  of  7 r0  (7To  — > ►  2j).  The  conditional  photon  energy  distribution  for  a  given  Ps 
energy  Eps  is  flat  and,  normalized,  is  given  by 


(16) 

where  the  photon  energy  e  is  between  the  Dopplcr-shift  limits 

7ps(1  -  Pps)rnec2  <  e  <  7Ps(l  +  /3Ps)mec2  (17) 

and  vanishes  outside  this  range.  Using  the  conditional  Ps  energy  distribution  for  a  positron 
of  energy  E  (Fps(Eps\E);  equation  12  in  §4)  and  integrating  over  the  Maxwellian  distribution 
of  the  positron,  the  photon  energy  distribution  for  H  is  then 


A(e)  oc  /  dEPsF,(e\Eps)  /  dE FPs{EPs\E)  VE e~E/kT 


6.8  eV+(e— mec2)2/e 


VEe-E/kTace(E)v(E) 

\E  —  6.8  eV)  (E  —  6.8  eV  +  4mec2) 


The  lower  limit  of  integration  is  the  minimum  energy  the  Ps  needs  to  produce  a  photon  of 
energy  e  (Stecker  1971).  For  He,  6.8  eV  is  replaced  by  17.8  eV.  In  Figure  8  we  show  the 
resulting  line  shapes  for  thermal  charge  exchange  in  H  (solid  curves)  and  He  (dashed  curves) 
for  various  temperatures. 


As  shown  in  §3,  the  energy  distributions  of  Ps  formed  in  flight  from  primary  positrons 
and  from  second-generation  positrons  resulting  from  ionization  of  primary  in-flight  Ps  are  also 
not  Maxwellian  (see  Figure  5).  From  our  Monte  Carlo  simulations  we  find  that  the  shapes 
of  the  associated  annihilation  lines  are  not  Gaussian.  The  solid  and  dashed  curves  in  Figure 
9  show  the  line  shapes  from  primary  and  second-generation  Ps,  respectively,  in  a  neutral 
medium  (He/H  =  0.1).  While  these  shapes  are  not  quite  Gaussian,  we  will  characterize  their 
widths  by  a  formal  full- width  at  half  maximum.  The  FWHM  of  the  primary  line  is  ~6.1  keV. 
(We  obtained  5.9  and  7.1  keV  FWHM  for  pure  H  and  He,  respectively.  The  value  of  5.9  is 
somewhat  less  than  the  value  of  6.4  keV  obtained  by  BRD  for  pure  H.)  If  the  medium  is  only 
slightly  ionized,  positrons  just  above  the  Ps  formation  threshold  are  rapidly  thermalized, 
slightly  modifying  the  Ps  formation  energy  distribution  (see  §3  and  dotted  histogram  of 
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Figure  5).  As  a  result,  there  are  fewer  photons  with  energies  very  near  511  keV,  flattening 
the  line  profile  as  shown  by  the  dotted  distribution  in  Figure  9.  At  ~6.9  keV,  the  FWHM  is 
slightly  larger  than  the  neutral  case.  The  FWHM  of  the  second-generation  line  is  ~2.6  keV 
(in  a  neutral  medium  with  He/H  =  0.1),  narrower  than  the  first-generation  line,  reflecting 
the  shift  to  lower  energies  of  the  Ps  energy  distribution  (see  Figure  5).  These  results  are 
only  weakly  dependent  on  density  and  temperature  (see  §3). 

To  obtain  the  shape  of  the  line  produced  by  direct  annihilation  of  thermal  positrons  with 
bound  electrons  it  is  necessary  to  know  the  momentum  distribution  of  the  atomic  electrons. 
We  first  treat  the  case  where  the  energies  of  the  positrons  are  near  zero.  In  this  case, 
the  shape  of  the  line  is  due  mostly  to  the  momentum  distribution  of  the  bound  electron. 
For  H,  this  distribution  was  calculated  by  Humberston  &  Wallace  (1972)  and  for  He  by 
Humberston  (1974)  and  Van  Reeth  et  al.  (1996).  These  calculations  showed  that  the  shape 
of  the  annihilation  line  is  nearly  Gaussian  with  a  FWHM  width  of  1.4  keV  for  H  and  2.5 
keV  for  He.  These  results  have  been  confirmed  by  laboratory  measurements  (Van  Reeth  et 
al.  1996). 

At  temperatures  greater  than  ~3xl04  K  for  H  and  ~lxl05  K  for  He,  the  Maxwellian 
momentum  distribution  of  the  thermal  positrons  begins  to  affect  the  shape  of  the  line.  At 
higher  temperatures,  the  effect  of  the  thermal  positron  dominates,  and  the  shape  of  the  line 
is  Gaussian.  The  shape  of  the  combined  line  for  all  temperatures  can  then  be  approximated 
by  a  Gaussian  with  a  FWHM  (r)  given  by 

r2  =  r2_  +  r2+  (19) 

where  re+  and  Te-  are  the  widths  for  the  thermal  positron  and  bound  electron  components, 
respectively,  with  re-  =  1.4  keV  for  H  and  2.5  keV  for  He  and 


derived  similarly  to  the  derivation  given  in  Appendix  A  for  the  shape  of  the  line  from  direct 
annihilation  of  positrons  and  free  electrons. 

We  summarize  the  various  shapes  and  widths  of  the  annihilation  line  from  the  five  modes 
of  annihilation. 

1.  Direct  annihilation  of  thermal  positrons  with  free  electrons.  The  shape  of 
the  line  is  Gaussian  with  the  FWHM  given  by  equation  (15).  The  annihilation  rate  by  this 
mode  is  shown  in  Figure  6. 

2.  Radiative  combination  of  thermal  positrons  with  free  electrons.  The  shape 
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of  the  line  from  the  xPs  decay  is  Gaussian  with  the  FWHM  given  by  equation  (15).  The  Ps 
formation  rate  by  this  mode  is  shown  in  Figure  6. 

3.  Positronium  formation  in  flight.  The  shape  of  the  line  from  the  xPs  decay  is  not 
Gaussian  (see  Figure  9)  and  has  a  FWHM  of  ~6.1  keV  for  Ps  from  primary  positrons  and 
~2.6  keV  for  Ps  from  second-generation  positrons.  The  fractions  of  primary  and  second- 
generation  positrons  forming  Ps  in  flight  are  shown  in  Figure  4. 

4.  Charge  exchange  of  thermal  positrons  with  atoms.  The  shape  of  the  line 

from  the  xPs  decay  is  not  Gaussian.  It  is  given  by  equation  (18)  and  is  shown  in  Figure  8. 
The  Ps  formation  rate  by  this  mode  is  shown  in  Figure  6. 

5.  Direct  Annihilation  of  Thermal  Positrons  with  Bound  Electrons.  The 

shape  is  approximately  Gaussian  with  the  FWHM  given  by  equation  (19).  The  annihilation 
rate  by  this  mode  is  shown  in  Figure  6. 


5.2.  Positronium  Continuum 


The  decay  of  3Ps  results  in  three  photons  which,  in  the  rest  frame  of  the  Ps,  are  dis¬ 
tributed  in  energy  according  to  (Ore  &  Powell  1949) 
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6.  Positronium  and  Thermal  Positron  Equilibrium  in  a  Plasma  and  the  Total 

Annihilation  Radiation  Spectrum 

As  shown  in  §5,  the  shape  of  the  2y  line  can  be  different  for  each  of  the  annihilation 
modes.  The  total  annihilation  radiation  spectrum  (both  2q  line  and  3q  continuum)  is  com¬ 
posed  of  contributions  from  each  of  the  modes.  In  order  to  construct  this  total  spectrum,  we 
need  to  know  the  relative  strengths  of  each  mode  for  a  given  combination  of  temperature, 
density  and  ionization  state  of  the  environment.  These  strengths  can  be  determined  rela¬ 
tively  easily  if  the  system  is  in  a  steady  state  condition.  It  will  be  in  steady  state  if  the  time 
scale  over  which  positrons  slow  down  and  annihilate  is  much  shorter  than  the  time  scale  over 
which  they  are  produced. 

The  time  history  for  positron  production  is  obtained  by  convolving  the  time  history 
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for  production  of  the  positron  emitters  with  their  exponential  decays.  The  decay  times 
of  the  most  important  positron  emitters  are  on  the  order  of  a  few  minutes  (Kozlovsky  et 
al.  1987).  The  time  required  for  positrons  to  slow  down  and  annihilate  can  be  estimated 
from  comparisons  of  calculated  positron  production  time  histories  with  measured  511  keV 
annihilation-line  time  histories  from  solar  flares.  Since  no  appreciable  differences  were  found 
(Murphy  et  al.  1984),  the  slowing-down  times  must  be  shorter  than  the  temporal  resolution 
of  the  instrument  (16  seconds  for  SMM)  which  is  shorter  than  the  positron  production  time 
scales.  The  assumption  of  steady  state  is  therefore  valid. 

In  §6.1  we  develop  a  system  of  conservation  equations  and  solve  it  for  the  numbers  of 
positrons  and  xPs  and  3Ps  atoms,  assuming  steady  state.  In  §6.2,  we  use  the  relative  numbers 
to  construct  the  total  annihilation  radiation  spectrum. 


6.1.  System  of  Continuity  Equations 


We  write  a  network  of  continuity  equations  for  the  conservation  of  the  number  densities 
of  thermal  positrons  and  Ps  formed  via  the  various  modes  and,  assuming  steady  state,  set 
all  of  their  partial  derivatives  with  respect  to  time  equal  to  zero.  The  2q  and  3q  production 
rates,  and  hence  the  emergent  spectrum,  can  then  be  obtained  from  the  resulting  number 
densities  which  all  scale  linearly  with  the  positron  injection  rate  per  unit  volume  Q+. 

The  network  of  reactions  was  displayed  graphically  in  Figure  1  and  particle  conservation 
is  expressed  through  the  following  system  of  continuity  equations: 
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Here  N+  is  the  number  of  thermal  positrons.  N{  and  N3  are  the  numbers  of  Ps  atoms  in  the 
singlet  and  triplet  states,  with  the  superscript  i  referring  to  the  Ps  formation  process:  Hce  - 
charge  exchange  with  H;  Hece  -  charge  exchange  with  He;  rc  -  radiative  combination  with 
free  electrons;  ifee  -  first-generation  in-flight  charge  exchange;  and  2 g  -  second-generation 
in-flight  charge  exchange;  each  of  which  has  its  own  Ps  velocity  distribution.  Rates  for  Ps 
formation  processes  (in  units  of  s^1)  are  denoted  Rl  where  the  superscript  i  has  the  same 
meaning  as  above.  These  rates  are  the  rates  calculated  above  (the  A’s,  see  §4)  multiplied 
by  the  relevant  target  densities.  Direct  annihilation  rates  are  also  denoted  Rl,  but  with  the 
superscript  i  referring  to  the  annihilation  channel:  Hdab  -  direct  annihilation  with  electrons 
bound  in  H;  Hedab  -  direct  annihilation  with  electrons  bound  in  He;  and  daf  -  direct 
annihilation  with  free  electrons.  Triplet  Ps  quenching  rates  are  denoted  R\^0  for  triplet- 
state  ionization,  R\_>i  for  triplet-singlet  spin  flip,  and  R\^x  =  Rs_> 0  +  where  again 

the  superscript  i  has  the  same  meaning  as  above. 

This  system  of  equations  can  then  be  solved  for  the  particle  numbers  .  First, 
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where  i  =  {Hce,  Hece ,  rc,  Hdab ,  Hedab,  daf}  and  j  =  {Hce,  Hece ,  rc}.  The  rest  of  the 
particle  numbers  can  be  expressed  in  terms  of  Q+  and  N+.  We  have 
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These  particle  numbers  depend  on  the  rates  R  which,  in  turn,  depend  on  the  density,  tem¬ 
perature  and  ionization  state  of  the  environment. 

From  the  particle  numbers  we  can  determine  the  rate  for  producing  2y  line  photons, 
Q27  =  2 N+  (i RHdab  +  RHedab  +  Rdaf^j  +  y  (N\fce  +  N^ce  +  N?ece  +  N[c  +  N ^ ,  (25) 

and  3y  continuum  photons, 


g37  =  1  (Nlfce  +  Af  ce  +  N*ece  +  iV3c  +  iVg9) .  (26) 

r3  V  / 

These  rates  are  directly  proportional  to  Q+  and  the  ratio  Q^/Qz-y  is  therefore  independent 
of  Q+. 


6.2.  Composite  Annihilation  Radiation  Spectrum 

The  composite  2y  annihilation  line  is  a  combination  of  the  line  shapes  from  each  of  the 
annihilation  modes  (see  §5.1).  The  relative  strength  of  each  mode  is  given  by  the  individual 
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terms  of  equation  (25)  as  calculated  for  the  temperature,  density  and  ionization  state  of  the 
environment.  The  composite  line  shape  can  be  very  different  from  the  shape  from  any  one 
of  the  modes. 

Let  Cr(e|eo,r)  denote  a  normalized  Gaussian  centered  at  eo  with  FWHM  denoted  by 
T.  Similarly,  let  H(e\eo,T)  denote  the  line  shape  resulting  from  thermal  charge  exchange 
and  J^ce(e|e0,r)  and  J29(e|e0,T)  denote  the  line  shapes  resulting  from  primary  and  second- 
generation  in-flight  charge  exchange,  respectively  (see  §5.1).  (We  note  that  these  latter 
shapes  can  be  far  from  Gaussian  and  the  formal  FWHM  may  not  adequately  reflect  this 
shape.)  Then  the  composite  2  7  line  shape  is  given  by 

^(e)  =  2N+[RHdabG(e\mec2,THdab) 

+RHedabG(e\meC2,THedab ) 


+RdafG(e\mec2,Tdaf) 


+N29J2°(e\mec2,T2g) 

+N?ceH{e\mec2,THce ) 

+N*eceH(e\mec2,  VHece) 

+Nr1cG{e\mec2,Trc)  . 

The  composite  3y  annihilation  continuum  is  composed  of  contributions  from  each  of  the 
four  Ps-prodncing  processes.  The  total  continuum  is  given  by 

^(<0  =  y /*'  [yyyy1^]  |Af-/“^“(e|£',ri/«)  +  Af32v2»(£|e',r29) 
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where  we  have  convolved  the  distribution  of  equation  (21)  with  the  line  shapes  for  each  mode 
to  approximate  the  effect  of  Doppler  broadening. 

7.  Applications 

We  will  apply  the  results  of  the  previous  sections  to  positron  annihilation  in  four  differ¬ 
ent  environments.  We  will  consider  a  completely  ionized  medium  (§7.1),  a  completely  neutral 
medium  (§7.2),  and  the  intermediate  case  of  a  partially- ionized  medium  (§7.3).  These  cases 
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will  provide  insight  on  how  the  line  shape  and  3Ps  continuum  depend  on  the  physical  condi¬ 
tions  of  the  medium.  This  insight  will  help  understand  the  results  when  we  then  consider  the 
fourth  environment:  a  specific  model  solar  atmosphere  (model  C  of  Vernazza  et  al.  (1981); 
§7.4)  which  provides  a  fixed  relation  between  temperature,  density  and  ionization  fractions. 

The  important  diagnostic  features  of  annihilation  radiation  are  the  Q?,1/Q2'1  ratio  and 
the  shape  of  the  line.  We  will  characterize  the  shape  of  the  line  with  its  full-width  at 
half  maximum  (FWHM).  However,  because  the  line  is  a  composite  of  the  many  modes 
of  annihilation,  its  shape  can  be  far  from  Gaussian,  and  we  caution  that  such  a  simple 
characterization  may  not  fully  represent  its  complexity.  In  particular,  when  the  broad  in¬ 
flight  component  is  present,  we  will  also  use  the  the  ratio  of  the  broad  component  yield  to  the 
sum  of  the  yields  of  the  narrow  components,  Rb/n,  as  an  additional  measure  of  the  line  shape. 
In  the  following  we  show  how  the  line  shape  and  Q^1/Q2l  depend  upon  the  conditions  of  the 
medium  where  the  annihilation  occurs.  The  medium  is  characterized  by  the  temperature 
T,  the  total  (neutral  +  ionized)  H  density  Wh,  and  the  ionization  fractions  XH+,  XHe+  and 
XHe++-  The  He  abundance  is  fixed  at  nue/nu  =  0.1  (Anders  &  Grevesse  1989). 


7.1.  Fully  Ionized  Medium 

We  set  XH+  =  XHe++  =  1,  XHe+  =  0  and  calculate  the  line  shape  and  Q^/Q-2^  for 
several  independent  choices  of  T  and  wr.  In  such  a  medium,  annihilation  can  occur  only  by 
direct  annihilation  with  free  electrons  and  Ps  formation  through  radiative  combination  with 
free  electrons.  Both  of  these  processes  produce  a  2y  line  with  a  Gaussian  shape  whose  width 
depends  only  on  temperature  and  is  given  by  equation  (15)  (see  §5.1). 

We  show  Q^/Q2^  as  a  function  of  T  for  several  densities  in  Figure  10.  The  decrease 
of  the  ratio  up  to  a  few  times  105  K  is  due  to  the  increasing  rate  of  positronium  quenching 
with  temperature  by  electrons  and  H+  (see  panel  (a)  of  Figure  7).  At  higher  temperatures, 
the  decrease  is  because  2y  annihilation  with  free  electrons  dominates  over  Ps  formation  via 
radiative  combination  (see  Figure  6). 

In  Figure  11,  Qz1/Q2'y  is  plotted  against  «h  for  several  temperatures.  For  each  temper¬ 
ature,  the  ratio  is  constant  for  densities  <1013.  This  is  because  at  these  densities  quenching 
is  not  important  and  both  Ps  formation  through  radiative  combination  with  free  electrons 
and  direct  annihilation  with  free  electrons  are  proportional  to  the  electron  density. 

The  profile  of  the  2y  line  and  3q  continuum  for  T  =  104  and  106  K  with  nH  =  1013 
cm-3  are  shown  in  panels  (a)  and  (b)  of  Figure  12,  respectively.  The  line  is  Gaussian  with 
width  given  by  equation  (15)  and  the  contribution  from  the  two  components  are  shown. 
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At  lower  temperature,  radiative  combination  dominates  while  at  higher  temperature,  direct 
annihilation  dominates. 

In  Figure  13  we  plot  the  Q37/Q27  ratio  versus  the  FWHM  line  width  for  several  densities. 
Because  both  of  the  processes  responsible  for  the  line  produce  Gaussian  line  shapes  with 
widths  depending  only  on  temperature  and  given  by  equation  (15),  the  widths  correspond 
directly  to  temperatures  and  the  corresponding  temperature  scale  is  shown  at  the  top  of  the 
figure. 


7.2.  Completely  Neutral  Medium 

In  a  completely  neutral  medium  (XH+  =  XHe+  =  Xne++  =  0  ),  annihilation  only  occurs 
via  Ps  (formed  by  in-flight  and  thermal  charge  exchange)  and  via  direct  annihilation  with 
bound  electrons.  At  the  low  temperatures  consistent  with  the  assumption  of  a  neutral 
medium,  the  contribution  from  thermal  charge  exchange  is  small  because  of  its  high  6.8  eV 
threshold  (the  high-energy  tail  of  the  positron  Maxwellian  is  too  weak).  The  Q^1/Q2^  ratio 
for  T  =  4000,  5000  and  6000  K  is  shown  as  a  function  of  density  by  the  dashed  curves  in 
Figure  11.  At  densities  greater  than  ~1013  cm-3,  the  ratio  falls  with  increasing  density  due 
to  quenching  of  3Ps  by  H.  At  densities  <1013  cm-3,  quenching  is  not  important  and  the  ratio 
is  constant.  At  4000  K,  for  example,  Q37/Q27  =  3.7.  Because  the  contribution  from  direct 
annihilation  with  bound  electrons  at  such  a  low  temperature  is  negligible,  this  value  can  be 
directly  derived.  From  Figure  4,  the  fraction  of  positrons  forming  positroninm  in  flight  is 
fi  =  0.95  for  Ah+  =  0  so  that  the  fraction  of  positrons  annihilating  directly  with  bound 
electrons  is  1  —  fy.  Therefore,  Q^/Qz-y  =  2.25/i/(0.5/i  +  2(1  —  fi ))  =  3.7  for  f\  =  0.95. 

The  composite  line  profile  (along  with  the  2y  line  components  and  the  3y  3Ps  continuum) 
is  shown  in  panel  (a)  of  Figure  14  for  T  =  5000  K  and  a  low  density  of  nn  =  109  cm-3  where 
quenching  is  not  important.  At  this  temperature,  the  2 7  annihilation  line  is  composed  mostly 
of  contributions  from  direct  annihilation  with  bound  electrons  and  a  strong  contribution  from 
charge  exchange  in  flight. 

At  higher  densities,  3Ps  quenching  modifies  the  line  shape  and  intensity.  Quenching 
of  3Ps  produced  by  thermalized  positrons  effectively  removes  three  photons  from  the  3 7 
continuum  and  adds  two  photons  to  the  narrow  2y  lines  (<2  keV  FWHM  at  the  temperatures 
consistent  with  a  neutral  assumption).  Quenching  of  primary  in-flight  3Ps  via  spin-flip  adds 
two  photons  to  the  broad  (~6.1  keV)  in-flight  2y  line.  Quenching  of  primary  in-flight  3Ps 
via  ionization  produces  second-generation  positrons  which  either  thermalize,  contributing 
to  the  narrow  2y  line  components,  or  form  Ps  again  in  flight.  This  second-generation  in- 


flight  Ps  either  annihilates  (to  produce  the  intermediate-width  ~2.6  keV  line  (see  §3.2)  and 
37  continuum)  or  is  again  quenched.  Because  most  of  the  Ps  atoms  formed  from  primary 
positrons  have  energies  greater  than  10  eV  (see  Figure  5),  quenching  via  ionization  is  favored 
over  spin-flip  due  to  their  respective  cross  sections  (see  Figures  3  and  panel  (d)  of  Figure  7). 
Therefore,  the  overall  effect  of  quenching  of  primary  in-flight  Ps  is  to  increase  the  yield  of  the 
narrower  lines  relative  to  the  broad  in-flight  line.  Panel  (b)  of  Figure  14  shows  the  calculated 
line  profile  for  T  =  5000  K  and  Uh  =  1017  cm-3.  Comparing  with  panel  (a)  for  109  cm-3,  the 
3y  continuum  is  now  absent,  the  narrow  components  from  thermal  positrons  annihilating 
via  charge-exchange  Ps  and  directly  with  bound  electrons  are  much  stronger,  the  intensity  of 
the  broad  in-flight  component  is  somewhat  increased  (due  to  spin-flip  quenching),  and  there 
is  now  a  contribution  from  second-generation  in-flight  charge  exchange. 

For  the  high-density  case  of  panel  (b),  the  line  is  dominated  by  the  narrow  components 
and  the  formal  FWHM  is  ~1.7  keV  (see  §5.1),  reflecting  only  these  narrow  components. 
For  the  low-density  case  of  panel  (a),  the  broad  in-flight  component  is  much  stronger,  yet 
the  formal  FWHM  is  still  only  ~2.6  keV,  intermediate  between  the  widths  of  the  narrow 
lines  and  the  broad  line  and  reflecting  neither.  These  examples  clearly  show  that,  except  for 
the  fully-ionized  case  where  the  broad  in-flight  component  is  never  present,  characterizing 
the  shape  of  the  annihilation  line  with  a  simple  measure  of  its  “width”  such  as  FWHM 
may  not  reveal  the  important  processes  or  the  true  shape  of  the  line.  Even  when  the  broad 
in-flight  component  is  relatively  weak,  as  in  the  high-density  neutral  case  of  panel  (b),  its 
presence  can  be  important  when  fitting  high  spectral  resolution  spectral  data  with  sufficiently 
good  statistical  significance.  When  the  in-flight  component  is  present,  the  line  is  essentially 
composed  a  broad  line  and  a  narrow  line.  In  this  case,  a  better  representation  of  the  shape 
of  the  line  is  Rb/n;  the  ratio  of  the  broad  component  yield  to  the  sum  of  the  yields  of  the 
narrow  components.  The  curves  in  Figure  15  show  Q^/Q2^  versus  Rb/n  for  T  =  4000,  5000 
and  6000  K  and  for  various  densities  as  noted  along  each  temperature  curve.  The  trend  of 
decreasing  and  Rb/n  with  increasing  density  along  each  curve  is  due  to  increased 

quenching:  3y  photons  are  moved  into  the  2y  line  which  decreases  arid  the  resulting 

build-up  of  the  narrow  2y  line  relative  to  the  broad  in-flight  line  decreases  Rb/n-  The  trend 
of  decreasing  arid  Rb/n  with  decreasing  temperature  for  a  given  density  is  due  to 

the  shift  from  Ps  formation  via  thermal  charge  exchange  to  direct  annihilation  with  bound 
electrons  as  temperature  decreases. 
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7.3.  Partially-Ionized  Medium 

Increasing  the  ionization  fraction  of  the  medium  decreases  Ps  formation  via  charge 
exchange  (by  both  in-flight  and  thermalized  positrons)  because  of  the  decreased  availability 
of  neutral  atoms.  In-flight  Ps  formation  is  further  reduced  because  the  larger  positron  energy 
loss  resulting  from  increased  electron  density  causes  more  positrons  to  thermalize  rather  than 
undergo  charge  exchange  in  flight.  This  is  reflected  in  a  decrease  of  both  f\  (the  fraction 
of  positrons  forming  Ps  in  flight)  and  Q^j Q^-  The  effect  on  the  line  profile  is  a  reduction 
of  the  broad  in-flight  component  and  a  corresponding  increase  of  the  narrow  component 
resulting  in  a  smaller  Rb/n-  Line  profiles  are  shown  in  Figure  16  for  several  values  of  XH+ 
(0.0,  0.01,  0.10  and  0.50)  for  T  =  5000  K  and  riu  =  1015  cm'3.  The  inset  table  gives  the 
associated  values  of  /i,  Qs1/Q2-y,  formal  line  width  and  Rb/n- 

We  have  compared  our  results  with  the  recent  positron  annihilation  calculations  of 
Churazov  et  al.  (2005)  for  a  low-density,  pure  H  medium  with  T  =  8000  K  and  XH+  =  0.1. 
For  the  line  shapes  we  find  good  agreement:  a  width  of  the  broad  line  from  in-flight  charge 
exchange  of  ~6  keV  (compared  with  ~5.9  keV  here)  and  a  width  of  the  narrow  line  from 
thermal  charge  exchange  of  ~  1.4  keV  (compared  with  ~1.3  keV  here).  However,  we  disagree 
on  fi,  the  fraction  of  positrons  that  form  Ps  in  flight.  While  Churazov  et  al.  (2005)  do  not 
provide  values  for  f\ ,  they  show  in  their  Figure  11  the  fraction  of  the  total  line  radiation 
contributed  by  the  broad  in-flight  line  which  is  directly  related.  Our  values  for  this  fraction 
are  consistently  higher  than  theirs;  for  example,  at  A"h+  =  0.1,  they  found  ~  0.31  for  this 
fraction  compared  with  ~0.52  from  our  calculations.  This  is  similar  to  the  difference  between 
our  calculations  for  f\  and  those  of  BRD  (see  §3.1). 


7.4.  Solar  Atmosphere  Model 

We  now  calculate  the  annihilation  line  shape  and  Qj,-y/Q2'y  throughout  the  fourth  envi¬ 
ronment:  the  quiet  solar  atmosphere  represented  by  the  non-LTE  calculations  of  Vernazza 
et  al.  (1981)  (model  C).  The  density,  temperature  and  ionization  fraction  vary  over  wide 
ranges  in  such  a  model  and  it  will  be  instructive  to  see  how  the  annihilation  radiation  varies 
in  response.  We  caution  that  atmospheric  conditions  during  a  solar  flare  may  be  quite  differ¬ 
ent  than  those  of  such  a  model.  Indeed,  as  we  will  see  in  §8,  observations  of  the  annihilation 
line  from  solar  flares  are  revealing  an  atmosphere  that  can  be  surprisingly  different.  We  also 
caution  that,  in  an  atmosphere  such  as  Vernazza  et  al.  (1981),  temperatures  and  ionization 
fractions  can  change  rapidly  over  very  narrow  ranges  of  height  (and  corresponding  gram- 
mage).  As  a  result,  at  any  given  point  there  may  not  be  sufficient  material  for  all  of  the 
positrons  to  remain  near  that  point  as  they  slow  down  and  annihilate.  Nevertheless,  the 
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calculation  will  demonstrate  the  range  of  and  line  shapes  possible  in  such  a  model 

atmosphere.  An  accurate  calculation  of  positron  annihilation  must  consider  the  location  of 
positron  emitter  production  and  transport  of  the  positrons,  with  the  result  that  the  line  may 
be  composed  of  contributions  from  a  range  of  conditions  in  the  solar  atmosphere. 

In  Figure  17  we  show  nH  and  ne  as  a  function  of  temperature  for  the  Vernazza  model 
using  the  tabulated  values  of  T,  nH  and  ne  (Vernazza  et  al.  1981).  (The  numbers  along  the 
H  curve  indicate  regions  that  will  be  referred  to  in  the  following  discussions.)  Vernazza  et  al. 
(1981)  do  not  specify  the  ionization  fractions  for  H  and  He  (i.e.,  XH+,  XHe+  and  XHe++).  We 
estimate  these  parameters  in  the  following  manner.  Up  to  the  temperature  where  ne  =  n h 
(T  =~  2.1  xlO4  K),  we  assume  that  He  is  totally  neutral;  i.e.,  the  free  electrons  are  due  solely 
to  ionized  H.  Above  this  temperature  we  assume  that  H  is  totally  ionized  and  the  excess  ne 
is  due  solely  to  singly-ionized  He  (He+).  Above  T  >  3.7  x  104  K  (at  which  ne/nu  =  1.1),  we 
assume  that  there  is  no  neutral  He  (He/H  is  assumed  to  be  10%).  We  thus  have  three  distinct 
temperature  regimes  where  the  ionization  fractions  have  the  values  (XH+,0,0),  (1,  VHe+,0) 
and  (1, 1  —  A7He++ ,  -VHe++ )  •  For  a  given  temperature,  then,  there  is  only  one  free  ionization 
parameter  which  is  set  by  the  ratio  ne/n h  tabulated  by  Vernazza  et  al.  (1981).  All  of  these 
parameters  can  be  related  to  the  height  in  the  atmosphere  using  Table  12  of  Vernazza  et  al. 
(1981). 

In  Figure  18,  Qz1/Q2l  is  plotted  as  a  function  of  temperature  (with  density  increasing 
along  the  curve  as  shown  by  the  arrow).  At  the  region  labeled  (1),  the  conditions  are  similar 
to  the  completely  neutral  case  discussed  above,  with  density  decreasing  rapidly  from  ~  1015 
to  ~  1013  cm-3.  As  quenching  decreases,  the  ratio  increases  until  it  reaches  the  maximum  of 
~4  (see  the  neutral  curve  of  Figure  11  for  T  =  6000  K).  At  region  (2),  H  is  partially  ionized 
and  He  is  neutral.  At  these  temperatures  and  ionization  states,  annihilation  is  dominated 
by  Ps  formation  (via  radiative  combination  and  both  thermal  and  in-flight  charge  exchange) 
and  the  ratio  reaches  its  maximum  possible  value  of  4.5  (i.e.,  (0.75x3)  /  (0.25x2))  since  the 
density  is  too  low  for  quenching  to  be  important.  Along  region  (3)  to  (4),  Q37/Q27  drops 
rapidly  due  to  H  becoming  fully  ionized  and  the  ratio  reaches  its  value  of  <3  for  fully  ionized 
as  shown  in  Figure  10  for  nn  —  1010  cm-3  and  T  ~  2  x  104  K.  In  region  (6),  there  are 
no  neutral  atoms  and  any  Ps  formation  is  due  to  radiative  combination  with  free  electrons. 
Q37/Q27  continues  to  decrease  due  to  increased  direct  annihilation  with  free  electrons  as 
temperature  increases  (see  Figure  6).  The  same  labeling  of  the  different  regions  is  employed 
in  Figure  19  where  is  plotted  versus  density.  We  see  that  for  such  an  atmosphere, 

a  given  value  of  Q^jQ^  can  occur  at  several  values  of  density. 

In  Figure  20  we  show  calculated  line  shapes  at  several  locations  in  the  atmosphere. 
In  panel  (a),  T  =  5000  K  and  riw  —  3^xl014  cm-3,  corresponding  to  region  (1)  of  Figure 
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17.  The  narrow  component  is  due  primarily  to  direct  annihilation  with  bound  electrons 
and  is  strong  relative  to  the  broad  component  due  to  Ps  formed  in  flight.  3Ps  quenching 
results  in  the  contribution  from  second-generation  in-flight  Ps.  In  panel  (b)  of  Figure  20, 
T  =  6000  K  and  %  —  2xl013  cm-3  and  the  broad  component  from  Ps  formed  in  flight 
reaches  its  relative  maximum  (the  narrow  peak  between  regions  (1)  and  (2)  of  Figure  18). 
The  narrow  component  now  has  a  significant  contribution  from  thermal  charge  exchange 
due  to  the  increasing  temperature  and  quenching  is  no  longer  significant.  In  panel  (c)  of 
Figure  20,  T  =  7000  K,  corresponding  to  %  —  1011  cm-3  and  region  (2)  of  Figure  18, 
and  the  broad  in-flight  charge  exchange  component  is  reduced  due  to  increased  ionization. 
The  narrow  component  is  due  to  thermal  charge  exchange  and,  because  of  the  increased 
ionization,  radiative  combination  and  direct  annihilation  with  free  electrons.  In  panel  (d) 
of  Figure  20,  T  =  104  K,  corresponding  to  %  —  5  x  1010  cm-3  (region  (3)  of  Figure  18). 
The  broad  component  is  negligible  and  the  narrow  component  is  due  primarily  to  thermal 
charge  exchange  with  the  remaining  neutral  atoms.  At  higher  temperatures  (lower  densities), 
thermal  charge  exchange  is  negligible  due  to  the  high  degree  of  ionization,  and  the  line  is 
formed  by  direct  annihilation  and  radiative  combination  with  free  electrons.  The  lines  are 
Gaussian  with  the  FWHM  given  by  equation  (15).  The  calculated  line  shapes  at  105  and 
106  K  are  shown  panels  (e  &  f)  of  Figure  20. 

In  Figures  21  and  22  we  plot  the  formal  FWHM  versus  density  and  temperature,  re¬ 
spectively  and  in  Figure  23  we  plot  the  ratio  versus  formal  FWHM.  In  all  three  of 

these  figures  we  also  plot  the  broad/narrow  ratio,  Rb/n-  When  Rb/n  is  greater  than  ~0.05, 
the  formal  FWHM  no  longer  reliably  reflects  the  shape  of  the  line  and  we  emphasize  this 
by  making  the  FWHM  curve  dotted  in  this  region.  Detailed  measurements  of  the  line  with 
high-resolution  detectors  such  as  RHESSI  can  reveal  the  true  shape  of  the  line  and  alter 
conclusions  based  on  FWHM  alone.  We  note  that  only  in  the  high-temperature/low-density 
regions  (regions  3-6  of  Figure  21)  does  the  formal  FWHM  ever  exceed  the  value  of  ~1.7  keV. 
We  also  see  that  the  relationships  are  not  unique,  so  that  measurement  of  just  one  quantity 
does  not  determine  where  in  the  atmosphere  the  annihilations  are  taking  place. 


8.  Application  to  RHESSI  Data 

RHESSI  (Lin  et  al.  2002)  made  the  first  high-resolution  spectral  observation  of  the 
solar  positron  annihilation  line  during  the  2002  July  23  solar  flare  (Share  et  al.  2003).  The 
line  was  subsequently  observed  in  both  the  2003  October  28  and  November  2  flares  (Share  et 
al.  2004).  When  Share  et  al.  (2004)  fit  the  annihilation  line  with  a  Gaussian,  the  measured 
width  during  all  three  flares  was  surprisingly  broad,  with  a  FWHM  of  4-8  keV.  The  October 
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28  flare  was  remarkable  in  that  during  the  first  10  minutes  of  the  observation  the  line  was 
broad  at  an  approximately  constant  value,  narrowed  dramatically  over  an  ~ 2-minute  interval 
to  ~1  keV,  and  then  remained  approximately  constant  at  this  value  during  the  remainder 
of  the  observation  (~12  minutes).  Except  for  the  Erst  four  minutes  of  the  October  28  flare, 
the  measured  Qz-y/Qz-y  ratio  in  all  of  the  flares  was  consistent  with  zero. 

In  the  previous  analyses,  the  data  were  compared  with  calculations  obtained  from  the 
same  computer  codes  used  in  the  current  paper,  but  the  cross  sections  for  the  various  interac¬ 
tions  were  not  yet  updated  and  the  process  of  in-flight  charge  exchange  of  second-generation 
positrons  was  not  included.  Since  those  analyses,  RHESSI  has  observed  the  annihilation 
line  in  at  least  one  other  flare,  that  of  2005  January  20.  The  analyses  of  the  data  from 
all  of  these  flares  continue  as  improvements  in  the  analysis  techniques  become  available.  In 
this  section  we  compare  the  calculations  with  data  but  only  as  an  example  of  what  can  be 
done  with  positron  annihilation  line  measurements  rather  than  trying  to  provide  definitive 
answers.  We  will  consider  RHESSI  observations  of  only  the  October  28  and  November  2 
flares  because  both  of  these  flares  were  more  intense  than  the  July  23  flare,  providing  data 
with  significantly-improved  statistical  quality. 

We  begin  with  the  November  2  flare.  In  their  analysis  of  this  flare,  Share  et  al.  (2004) 
obtained  an  acceptable  fit  to  the  annihilation  line  spectrum  using  a  4.8  ±  0.5  keV  (FWHM) 
Gaussian.  Using  the  new  calculation,  we  hnd  that  a  Vernazza  quiet-Sun  atmosphere  at 
T  ~  5000  K  also  provides  an  acceptable  fit  (13%  probability  based  on  y2)  to  the  500  -  520 
keV  spectrum  shown  in  panel  (a)  of  Figure  24;  the  difierent  components  of  the  model  are 
plotted  in  panel  (a)  of  Figure  20.  However,  we  do  not  believe  we  can  confidently  conclude  that 
the  annihilation  occurred  in  this  environment.  The  0.01  gm  cm-2  overlying  column  depth 
and  ~  3  x  1014  cm-3  density  at  this  height  in  the  quiet-Sun  atmosphere  are  only  marginally 
sufficient  for  production  of  the  radioactive  positron  emitters  by  accelerated  particles  with 
a  typical  power  law  index  of  ~  4.  The  spectrum  for  this  flare  may  be  considerably  harder 
(Share  et  al.  2004),  implying  a  large  fraction  of  positrons  from  the  decay  of  pions  whose  high- 
energy  production  reactions  require  an  even  greater  column  depth.  Under  these  conditions 
a  significant  fraction  of  the  annihilation  line  could  be  Compton  scattered,  with  the  scattered 
continuum  contributing  to  the  apparent  Q^1/Q2'y  ratio.  The  measured  99%  upper  limit  on 
the  ratio  (1.3)  is  only  marginally  consistent  with  the  calculated  ratio  (1.1)  even  without  a 
contribution  from  Compton  scattering. 

We  plot  the  line  spectrum  accumulated  during  the  first  2-minutes  of  the  October  28  flare 
in  panel  (b)  of  Figure  24.  It  is  clearly  too  broad  to  be  consistent  with  the  calculated  shape 
produced  by  annihilation  in  a  quiet-Sun  atmosphere  at  5000  K.  The  only  way  to  obtain  a 
broader  line  at  temperatures  <  104  K  is  to  increase  the  relative  contribution  from  in-flight 
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charge  exchange.  The  maximum  contribution  in  the  Vernazza  atmosphere  occurs  at  a  tem¬ 
perature  of  6000  K  (Rb/n  —  3.9,  see  Figure  22  and  panel  (b)  of  Figure  20).  The  calculated 
shape  for  this  environment  (solid  curve)  provides  a  marginally-acceptable  fit  (~1.2%  confi¬ 
dence).  Other  considerations  exclude  this  origin,  however:  the  6  x  10-4  gm  cm-2  overlying 
column  depth  and  ~  2  x  1013  cm-3  density  are  much  too  small  for  sufficient  nuclear  reactions 
and  the  calculated  Q3l/Q2l  ratio  (3.5)  is  inconsistent  with  the  measured  ratio  (1.8  ±  0.4). 
We  also  fit  the  data  with  line  shapes  derived  assuming  the  FI  flare  atmosphere  of  Machado 
et  al.  (1980)  but  no  acceptable  fit  to  the  line  could  be  obtained  since  Rb/n  never  exceeds 
1.5  in  the  cool  regions  of  such  an  atmosphere.  The  alternative  scenario  is  annihilation  in 
a  high-temperature  region  where  the  line  is  a  thermally-broadened  Gaussian.  Panel  (b)  of 
Figure  24  shows  this  Gaussian  fit  (dotted  curve;  6.7  ±  1.5  keV  FWHM).  From  equation 
(15),  this  width  corresponds  to  T  =  3  —  4  x  105  K.  At  this  temperature,  H  and  He  would  be 
ionized  and  Figures  10,  11  and  13  can  be  used.  Figure  13  shows  that  the  densities  required 
for  adequate  positron  production  (>  1015  cm-3)  correspond  to  Q3l/Q2l  ratios  <0.2  which 
are  significantly  less  than  the  measured  value  (1.8  ±  0.4).  However,  such  small  values  of 
Q3l/Q2l  may  be  allowed  since  much  of  the  continuum  below  the  511  keV  line  may  be  due 
to  Compton  scattering  of  the  line  at  the  large  depths  where  pion-decay  positrons  annihilate. 
We  base  this  argument  on  the  hardness  of  the  accelerated  particle  spectrum  and  detection 
of  the  characteristic  pion-decay  7-ray  spectrum  observed  by  CORONAS-F  (Kuznetsov  et 
al.  2005).  We  therefore  confirm  the  Share  et  al.  (2004)  conclusion  that  positron  annihi¬ 
lation  during  the  early  part  of  the  October  28  flare  occurred  in  a  high  temperature  region 
( T  =  3  —  4  x  105  K)  at  densities  >  1015  cm-3. 

In  panel  (c)  of  Figure  24  we  show  that  an  ~1  keV  FWHM  Gaussian  (dotted  curve) 
provides  an  excellent  fit  to  the  narrow  line  observed  during  the  last  10  minutes  of  the  October 
28  flare  (Share  et  al.  2004)  when  the  measured  99%  upper  limit  for  the  Q3l/Q2l  ratio  is  0.3. 
Such  a  narrow  width  requires  both  that  the  temperature  be  less  than  104  K  and  that  the 
broad/narrow  ratio  Rb/n  be  sufficiently  small.  We  find  that  there  is  no  temperature  in  either 
the  Vernazza  atmosphere  or  a  neutral  medium  at  which  this  can  be  achieved  and  also  satisfy 
the  Q3^/Q2l  constraint  (see  Figures  15  and  18  &  22).  An  acceptable  fit  can  be  achieved, 
however,  by  increasing  the  ionization  fraction  which  reduces  Rb/n  (see  Figure  16).  Assuming 
T  =  5000  K  and  a  density  of  1015  cm-3  (which  satisfies  both  the  Q^/Q^  constraint  and 
the  nuclear  production  requirement),  we  find  that  an  ionization  fraction  of  at  least  20%  is 
required  to  provide  an  acceptable  fit  (>1%  confidence)  to  the  spectrum.  The  solid  curve  in 
panel  (c)  of  Figure  24  shows  this  fit. 

In  summary,  we  find  that  positron  annihilation  in  solar  flares  occurs  in  a  wide  range  of 
conditions.  The  line  shape  observed  during  the  November  2  flare  was  adequately  fit  with 
the  quiet  solar  atmosphere  of  Vernazza  et  al.  (1981)  at  T  =  5000  K.  When  the  line  was 
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broad  during  the  October  28  flare,  the  data  suggest  that  the  environment  was  ionized  with 
T  =  3  — 4  x  105  K  at  a  density  >  1015  cm-3.  When  the  line  was  narrow,  the  environment  was 
relatively  cool  at  ~5000  K  yet  the  H  ionization  fraction  had  to  have  been  at  least  20%.  We 
emphasize  again  that  these  results  are  meant  only  as  an  example  of  the  application  of  the 
calculations  to  data.  More  detailed  analyses  of  these  flares  will  be  presented  in  a  subsequent 
paper. 


9.  Summary 

The  formation  of  the  positron  annihilation  line  is  a  complex  interplay  of  competing  pro¬ 
cesses.  In  this  paper  we  explored  this  process,  identifying  the  various  annihilation  modes, 
determining  the  nature  of  the  contribution  of  each  mode,  and  establishing  how  the  phys¬ 
ical  conditions  of  the  environment  determine  their  relative  contributions.  We  thoroughly 
updated  the  cross  sections  for  all  of  the  relevant  processes  with  recent  laboratory  measure¬ 
ments  and  theoretical  calculations,  providing  reliable  cross  sections  for  all  of  the  important 
interactions.  We  calculated  the  relative  rates  for  each  of  the  modes,  their  associated  line 
shapes,  and  how  the  individual  contributions  combine  to  produce  the  final  composite  line 
shape.  We  considered  neutral,  fully  ionized,  and  partially  ionized  mediums,  and  the  quiet 
sun  atmosphere  of  Vernazza  et  al.  (1981). 

In  a  fully  ionized  medium  with  T  >  104  K,  Qz1jQ2^  is  less  than  ~3  and  decreases  with 
increasing  density  (due  to  collisional  quenching  of  3Ps)  and  increasing  temperature  (due  to 
direct  annihilation  dominating  over  Ps  formation  via  radiative  combination).  The  width  of 
the  line  depends  directly  on  temperature  as  given  by  equation  (15).  In  a  neutral  medium, 
Q37/Q27  reaches  its  maximum  value  at  low  densities  where  quenching  is  not  important  (for 
T  =  5000  K,  Q37/ Q27  —  3.6).  The  line  is  composed  of  a  narrow  component  (~  1.5  keV) 
due  to  direct  annihilation  with  bound  electrons  and  a  broad  component  (~  6.1  keV)  due  to 
charge  exchange  in  flight  with  neutral  H.  The  latter  can  contribute  significantly  to  the  line 
shape  when  the  density  is  low  and  quenching  is  not  important.  The  effect  of  increasing  the 
ionization  fraction  is  to  decrease  Ps  formation  via  charge  exchange  (by  both  in-flight  and 
therm ali zed  positrons)  due  to  the  decreasing  availability  of  neutral  atoms  which  reduces  the 
broad  component  relative  to  the  narrow  components. 

For  the  model  C  solar  atmosphere  of  Vernazza  et  al.  (1981),  we  calculated  the  annihi¬ 
lation  line  shape  and  the  Qz1IQ,2^  ratio  throughout.  Such  a  model  provides  a  wide  range  of 
physical  conditions  and  we  identified  how  these  conditions  affect  the  line  shape  and  Qz1/Q2~r 
With  such  an  understanding,  analyses  of  annihilation  line  observations  from  solar  flares  can 
provide  reliable  constraints  on  the  physical  conditions  where  the  positrons  annihilate. 
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We  compared  the  calculated  line  shapes  and  Q^IQi^  ratios  with  high  spectral  resolution 
data  recently  obtained  with  RHESSI  from  the  2003  October  28  and  November  2  flares.  We 
found  that  positron  annihilation  in  solar  flares  can  occur  in  a  wide  range  of  conditions.  The 
annihilation  environment  of  the  November  2  flare  could  have  been  at  a  temperature  of  5000 
K  in  the  quiet  solar  atmosphere  of  Vernazza  et  ah  (1981).  When  the  line  was  broad  during 
the  October  28  flare,  the  data  suggest  that  the  environment  was  ionized  with  T  =  3  — 4  x  105 
K  and  density  >  1014  cm-3.  When  the  line  was  narrow,  the  environment  was  relatively  cool 
at  ~5000  K  yet  the  H  ionization  fraction  had  to  have  been  at  least  20%. 


We  thank  H.  R.  J.  Walters  for  valuable  discussions  and  for  providing  unpublished  cross 
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A.  Direct  Annihilation  Line  Spectrum 


Here  we  show  that  the  line  resulting  from  the  direct  annihilation  of  nonrelativistic 
thermal  electrons  and  positrons  has  a  gaussian  distribution  with  width  given  by 


FWHM  =  mCyin2^  =l.lkeVy^.  (Al) 

In  a  free  annihilation  the  emitted  photon  energies  in  the  lab  system  are 

e=  (l+  E+  +  f^')  (1  ±  f}cmH*)mec 2  (A2) 

V  mecz  J 

where  E+  and  E _  are  the  kinetic  energies  of  the  positron  and  electron,  u  =  c/3cm  is  the 
magnitude  of  the  velocity  of  the  center  of  mass,  and  (i*  is  the  cosine  of  the  angle  between 
one  of  the  emitted  photons  and  the  electron  in  the  CM  system.  Again,  all  starred  quantities 
refer  to  the  CM  system.  At  nonrelativistic  energies  (E±  <C  mec2)  this  expression  reduces  to 

e  ~  (1  ±  / 3cmn*)mec 2.  (A3) 


The  cross  section  differential  in  the  cosine  of  the  polar  angle  fi*  is  (e.g.  Heitler  1954) 


da* 

dfi* 


7T  r-n 


4  7 *2(3* 


7*2  +  (7*2  _  1)(2  _  p*2)  2(7*2  -  1)2(1  -  fi*2)2 


7* 2  —  (7*2  —  l)^*2]^ 


(A4) 


7*2  —  (7*2  —  l)/i*2 
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In  the  nonrelativistic  limit  the  photons  are  emitted  isotropically  in  the  CM  system  because 
this  cross  section  reduces  to 


da* 


TTTqC 


dp*  2v  ’  ^ 

which  is  independent  of  angle  and  only  depends  on  the  relative  velocity  v  =  |v+  —  v_|.  This 
holds  even  with  the  inclusion  of  the  Coulomb  correction  factor  (2nac/v)/(l  —  e~27rac/v)  in 
which  case  the  cross  section  becomes 

da*  7rrg  2ir a(c/v)2 


dp* 

The  annihilation  spectrum  is  given  by 


2  ^  _  g— 2'Kac/v  ' 


(A6) 


dN~ 


dt  de  d3x 


(£)  =  2n+n-(^f)S  fd3v*d3v- 


exp 


mev+ 
'  2kT 


exp 


rnev2_ 
’  2kT 


da*  ,  . 

v- — (v) 


dfi* 


dji* 


de 


Equation  (A3)  gives 


dp* 


de 


rripcu 


(AT) 


(AS) 


Equation  (A3)  also  implies  that  to  obtain  a  photon  of  energy  e  the  minimum  CM  velocity  is 


u0(e)  =  c 


1 rripC * 


(A9) 


Under  the  coordinate  transformation  {v+,v_}  — >  {u,  v},  where  u  is  the  velocity  of  the 
center  of  mass,  the  annihilation  spectrum  becomes 


dN~ 


dt  de  d3x 


.e  = 


2n+n-  (2 I  Mvexp(- 


meu 

~kT 


2n+n_(47r)2  /  me  \3 
mec  \2nkT ) 


Co(e) 


dvv2  exp 
duu  exp 


exp 


mev 

'JkT 


mev2  \  da* 

w) 

da*  , 


dp* 


mecu 

(A10) 


mPu 


kT 


The  only  dependence  on  e  is  through  «o(e)  in  Uie  integration  over  u.  The  integration  over  v 
simply  yields  an  overall  constant,  hence  we  can  write 


1  gL/V7 
A"7  de 


°c  / 

J  uo(e) 
oc  exp 

=  exp 


duu  exp 

meul 

kT 

mec 2 


meu 

~kT 


(All) 


kT  V  mec 2 
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Finally,  from  (All)  it  is  a  simple  matter  to  show  that  the  full  width  at  half  maximum  is 
given  by  equation  (Al). 
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Table  1.  Interactions  Involving  Positrons  and  Positronium 


Energy  Redistribution  Positronium  Formation  Positronium  Quenching  Annihilation 


e+  H  ->  e+  H* 
e+  H  ->  e+e“  H+ 
e+  He  — >  e+  He* 
e+  He  -»•  e+e~  He+ 


e+e  — ■»  7  Ps 
e+  H  ->  H+  Ps 
e+  He  ->•  He+  Ps 


e  3Ps  - 

-y  e  e  e+ 

H  3Ps  - 

►  H  e“e+ 

1 

CO 

Oh 

CO 

1 

03 

->  e“  3Ps 

H  3Ps  - 

►  H  3Ps 

H+  Ps  - 

4  H  e+ 

He+  Ps  - 

-*•  He  e+ 

e+e  — >  2y 
e+  H  ->  H+  27 
e+  He  — >  He+  27 
xPs  -»•  27 
3Ps  — >•  37 
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Fig.  1  —  Schematic  representation  of  the  processes  by  which  positrons  annihilate. 

Fig.  2  —  Cross  sections  for  ionization  (ion),  excitation  (ex)  and  charge  exchange  (ce)  in¬ 
teractions  of  positrons  with  H  (solid)  and  He  (dashed).  Also  shown  is  the  cross  section  for 
direct  annihilation  of  positrons  with  free  electrons  (daf;  dotted),  multiplied  by  105. 

Fig.  3  —  Cross  sections  for  quenching  interactions  of  Ps.  If  the  interaction  is  with  an  atom 
or  ion  then  the  abscissa  gives  the  Ps  energy.  If  the  interaction  is  with  an  electron  then  the 
abscissa  gives  the  energy  of  the  electron  in  the  rest  frame  of  the  Ps. 

Fig.  4  —  The  fraction  of  positrons  which  form  Ps  in  flight  as  a  function  of  the  ionization 
fraction  of  H  (Ah+)  for  ne  —  5  x  1013  cm-3  and  T  =  11,600  K.  The  various  curves  are 
discussed  in  the  text. 

Fig.  5  —  The  Ps  energy  distribution  resulting  from  formation  in  flight  by  primary  positrons 
calculated  by  the  Monte  Carlo  simulation  for  AH+  =  0  (solid)  and  XH+  =  0.05  (dotted). 
The  dashed  histogram  show  the  distribution  for  secondary  positrons  resulting  from  ionization 
quenching  of  primary  3Ps  (AH+  =0). 

Fig.  6  —  Thermally  averaged  rates  per  unit  target  density  for  direct  annihilation  (daf)  and 
radiative  combination  (rc)  with  free  electrons,  charge  exchange  (ce)  with  H  and  He  and  direct 
annihilation  with  bound  electrons  (dab)  of  H  and  He.  The  curve  for  direct  annihilation  with 
bound  electrons  of  H  ends  at  T  ~  104  K  since  above  this  temperature,  H  is  considered  totally 
ionized. 

Fig.  7  —  Thermally  averaged  quenching  rates  calculated  from  the  cross  sections  of  Fig.  3. 
The  Ps  has  the  velocity  distribution  resulting  from  Ps  formation  through  thermal  radiative 
combination  in  panel  (a),  thermal  charge  exchange  with  H  in  panel  (b),  thermal  charge 
exchange  with  He  in  panel  (c),  and  charge  exchange  in  flight  in  panel  (d).  We  note  that  at 
the  higher  temperatures,  neutral  H  or  He  may  not  be  present  fo  charge  exchange  to  occur. 

Fig.  8  —  Line  profiles  for  the  2 7  annihilation  of  3Ps  formed  through  charge  exchange  of 
thermal  positrons  with  H  and  He  for  several  temperatures. 

Fig.  9  —  Line  profile  for  the  2y  annihilation  of  4Ps  formed  through  charge  exchange  in 
flight  using  the  Ps  energy  distributions  from  primary  positrons  shown  in  Figure  5:  XH+  =  0 
(solid)  and  XH+  =  0.05  (dotted).  The  dashed  histogram  shows  the  line  profile  using  the  Ps 
distribution  from  secondary  positrons  (A"H+  =  0). 

Fig.  10  —  Ratio  of  the  number  of  3y  continuum  photons  to  2y  line  photons  ( Q^IQ2^ )  as  a 
function  of  temperature  for  various  densities  in  a  fully  ionized  medium. 
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Fig.  11  —  Ratio  of  the  number  of  37  continuum  photons  to  27  line  photons  (Q^/Q 2 7)  as  a 
function  of  density  for  various  temperatures  in  a  fully  ionized  medium.  The  dashed  curves 
are  for  a  neutral  medium  at  T  =  4000,  5000  and  6000  K. 

Fig.  12  —  Calculated  line  profiles  for  fully  ionized  medium  at  T  =  104  and  106  K  with 
n h  =  1013  cm”3.  The  components  contributing  to  the  spectrum  are  indicated:  radiative 
combination  with  free  electrons  (rc),  direct  annihilation  with  free  electrons  (daf)  and  the  3y 
3Ps  continuum  (3y,  dotted  curve).  The  normalization  is  1  initial  positron. 

Fig.  13  —  Ratio  of  the  number  of  3y  continuum  photons  to  2y  line  photons  (Q37/Q27)  as 
a  function  of  line  width  for  various  densities  in  a  fully  ionized  medium.  Because  in  a  fully 
ionized  medium  the  width  is  directly  related  to  the  temperature  through  equation  (15),  we 
show  the  corresponding  temperature  scale  along  the  top  of  the  figure. 

Fig.  14  —  Calculated  line  profiles  for  a  neutral  medium  at  riu  =  109  and  101'  cm”3  with 
T  =  5000  K.  The  components  contributing  to  the  spectrum  are  indicated:  in-flight  charge 
exchange  (ifee),  direct  annihilation  with  bound  H  electrons  (dab,H),  direct  annihilation  with 
bound  He  electrons  (dab, He),  thermal  charge  exchange  with  H  (ce,H),  second-generation 
in-flight  charge  exchange  (2g),  and  3y  3Ps  continuum  (3j,  dotted  curve).  The  normalization 
is  1  initial  positron. 

Fig.  15  —  Ratio  of  the  number  of  3j  continuum  photons  to  2y  line  photons  (Q37/Q27)  in 
a  neutral  medium  as  a  function  of  the  broad/narrow  line  ratio  Rb/n  for  T  =  4000,  5000  & 
6000  K  and  various  densities  as  indicated  along  each  curve 

Fig.  16  —  Calculated  line  profiles  for  a  partially-ionized  medium  with  T  =  5000  K,  nH  = 
1015  cm”3,  and  XH+  =  0.0,  0.01,  0.10  and  0.50.  The  inset  table  gives  the  associated  values 
°f  /1,  Q37/Q27,  formal  line  width  (FWHM),  and  the  broad/narrow  line  ratio  Rb/n-  The 
normalization  is  to  1  initial  positron. 

Fig.  17  —  The  total  H  (H°  +  H+)  and  electron  densities  versus  temperature  in  the  solar 
atmosphere  according  to  model  C  of  Vernazza  et  al.  (1981).  The  numbers  along  the  curve 
correspond  to  regions  discussed  in  the  text. 

Fig.  18  —  Ratio  of  the  number  of  3y  continuum  photons  to  2y  line  photons  (Q37/Q27)  as 
a  function  of  temperature  with  the  density  and  ionization  fractions  constrained  by  model  C 
of  Vernazza  et  al.  (1981).  The  numbers  along  the  curve  correspond  to  the  regions  noted  in 
Figure  17. 

Fig.  19  —  Ratio  of  the  number  of  3y  continuum  photons  to  2y  line  photons  (Q37/Q27)  as 
a  function  of  density  with  the  temperature  and  ionization  fractions  constrained  by  model  C 
of  Vernazza  et  al.  (1981).  The  numbers  along  the  curve  correspond  to  the  regions  noted  in 
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Figure  17. 

Fig.  20  —  Line  profiles  for  various  temperatures  with  the  density  and  ionization  fractions 
constrained  by  model  C  of  Vernazza  et  al.  (1981).  The  annihilation  modes  are  indicated: 
in-flight  charge  exchange  (ifce),  thermal  charge  exchange  with  H  (ce,H),  thermal  charge 
exchange  with  He  (ce,He),  radiative  combination  with  free  electrons  (rc),  direct  annihilation 
with  bound  H  electrons  (dab,H),  direct  annihilation  with  bound  He  electrons  (dab, He),  direct 
annihilation  with  free  electrons  (daf),  second-generation  in-flight  charge  exchange  (2g),  and 
3y  3Ps  continuum  (3y,  dotted  curve).  The  normalization  is  to  1  initial  positron. 

Fig.  21  —  Line  FWHM  and  the  broad/narrow  line  ratio  Rb/n  as  functions  of  density  with 
the  temperature  and  ionization  fractions  constrained  by  model  C  of  Vernazza  et  al.  (1981). 
The  numbers  along  the  curve  correspond  to  the  regions  noted  in  Figure  17.  FWHM  is  on 
the  left  axis  and  Rb/n  is  on  the  right. 

Fig.  22  —  Line  FWHM  and  the  broad/narrow  line  ratio  Rb/n  as  functions  of  temperature 
with  the  density  and  ionization  fractions  constrained  by  model  C  of  Vernazza  et  al.  (1981). 
The  numbers  along  the  curve  correspond  to  the  regions  noted  in  Figure  17.  FWHM  is  on 
the  left  axis  and  Rb/n  is  on  the  right. 

Fig.  23  —  Ratio  of  the  number  of  3y  continuum  photons  to  2y  line  photons  as 

a  function  of  the  line  FWHM  and  the  broad/narrow  line  ratio  Rb/n  with  the  density  and 
ionization  fractions  constrained  by  model  C  of  Vernazza  et  al.  (1981).  The  numbers  along 
the  curve  correspond  to  the  regions  noted  in  Figure  17.  FWHM  is  on  the  bottom  axis  and 
Rb/n  is  on  the  top. 

Fig.  24  —  Observed  count  spectra  and  fits  for  the  2003  November  2  and  2003  October 
28  flares.  The  November  2  spectrum  is  shown  in  (panel  (a),  along  with  the  T  =  5000  K 
Vernazza  fit.  The  October  28  broad  line  is  shown  in  panel  (b),  overplotted  with  the  T  = 
6000  K  Vernazza  (solid  curve)  and  Gaussian  (dotted  curve)  fits.  The  October  28  narrow 
line  is  shown  in  panel  (c),  over  plotted  with  the  20%  ionized  (solid  curve)  and  the  Gaussian 
(dotted  curve)  fits. 
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